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Abstract 

This is the second in a series of papers in which we derive a T-expansion for the two- 
dimensional non-local Ginzburg-Landau energy with Coulomb repulsion known as the 
Ohta-Kawasaki model in connection with diblock copolymer systems. In this model, 
two phases appear, which interact via a nonlocal Coulomb type energy. Here we focus 
on the sharp interface version of this energy in the regime where one of the phases 
has very small volume fraction, thus creating small "droplets" of the minority phase 
in a "sea" of the majority phase. In our previous paper, we computed the T-limit of 
the leading order energy, which yields the averaged behavior for almost minimizers, 
namely that the density of droplets should be uniform. Here we go to the next order 
and derive a next order T-limit energy, which is exactly the Coulombian renormalized 
energy obtained by Sandier and Serfaty as a limiting interaction energy for vortices 
in the magnetic Ginzburg-Landau model. The derivation is based on the abstract 
scheme of Sandier-Serfaty that serves to obtain lower bounds for 2-scale energies and 
express them through some probabilities on patterns via the multiparameter ergodic 
theorem. Without thus appealing to the Euler-Lagrange equation, we establish for 
all configurations which have "almost minimal energy" the asymptotic roundness and 
radius of the droplets, and the fact that they asymptotically shrink to points whose 
arrangement minimizes the renormalized energy in some averaged sense. Via a kind of 
Inequivalence, the obtained results also yield an expansion of the minimal energy for 
the original Ohta-Kawasaki energy. This leads to expecting to see triangular lattices 
of droplets as energy minimizers. 



1 Introduction 

This is our second paper devoted to the T-convergence study of the two-dimensional Ohta- 
Kawasaki energy functional [28] in two space dimensions in the regime near the onset of 
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non-trivial minimizers. The energy functional has the following form: 



£[u] = I ( e -\Vu\ 2 + V(u)\dx + ^ I f ( U (x)-u)G (x,y)(u(y)-u)dxdy, (1.1) 

where f2 is the domain occupied by the material, u : VL — > R is the scalar order parameter, 
V(u) is a symmetric double-well potential with minima at u = ±1, such as the usual 
Ginzburg-Landau potential V(u) = ^(1 — u 2 ) 2 (for simplicity, the overall coefficient in V 
is chosen to make the associated surface tension constant to be equal to e, i.e., we have 
f \ sj2V{u) du = 1), e > is a parameter characterizing interfacial thickness, u E (— 1, 1) 
is the background charge density, and Go is the Neumann Green's function of the Laplacian, 
i.e., Go solves 

-AG (x,y) = S(x-y)--^-, G (x,y) dx = 0, (1.2) 

where A is the Laplacian in x and 5(x) is the Dirac delta-function, with Neumann boundary 
conditions. Note that u is also assumed to satisfy the "charge neutrality" condition 
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udx = u. (1-3) 
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For a discussion of the motivation and the main quantitative features of this model, see 
our first paper [19], as well as [25j[26]. For specific applications to physical systems, we 
refer the reader to PEIIIB1I221I231I251I23I2H10D] • 

In our first paper [19], we established the leading order term in the T-expansion of the 
energy in (jl.lh in the scaling regime corresponding to the threshold between trivial and 
non-trivial minimizers. More precisely, we studied the behavior of the energy as e — > 
when 

u £ := -l + e 2/3 |lne| 1/3 5, (1.4) 

for some fixed 5 > and when f2 is a flat two-dimensional torus of side length £, i.e., when 
fi = T? = [0, i) 2 , with periodic boundary conditions. As follows from [191 Theorem 2] and 
the arguments in the proof of [191 Theorem 3], in this regime minimizers of E consist of 
many small "droplets" (regions where u > 0) and their number blows up as e — > 0. We 
showed that, after a suitable rescaling the energy functional in (ll.ip T-converges in the 
sense of convergence of the (suitably normalized) droplet densities, to the limit functional 
E°[n] defined for all densities fi G M(Tj) n H~ l (T 2 ) by: 

E°[fA = Pj + ( 3 2/3 - %) [ du + 2 // G(x- y)dfJL(x)diM(y), (1.5) 
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where G(x) is the screened Green's function of the Laplacian, i.e., it solves the periodic 
problem for the equation 

-AG + k 2 G = 5(x) in Tj, (1.6) 

and k = l/y/V"(l) = §. Here we noted that the double integral in (|1.5|) is well defined 
in the sense fJ* T 2 xT 2 G(x — y)dfi(x)dfi(y) := L 2 vdfi, where the latter is interpreted as the 
Hahn-Banach extension of the corresponding linear functional defined by the integral on 
smooth test functions (see also [Ml Sec. 7.3.1] and [9] for further discussion). Indeed, 
v := G * dfj, is the convolution understood distributionally, i.e., (G * d/i, /} := (G * /, dfj) = 
fj 2 (^Jj2 G(x — y)f(y)dy S j dfi(x) for every / G C°°(T|) and, hence, by elliptic regularity 
IMI/f^Tf) — CII/lliJ-^T 2 ) for some C > 0, so v G H l (Tf). 
In particular, for 5 > 5 C , where 

S c := (1.7) 

the limit energy i? [/i] is minimized by d/j,(x) = fidx, where 

fi = l(S-S c ) and E [J1} = ^(25-6 C ). (1.8) 

When 5 < 5 C , the limit energy is minimized by [i = 0, with E°[0] = 5 2 /(2k 2 ). The value of 
5 = S c thus serves as the threshold separating the trivial and the non-trivial minimizers of 
the energy in (jl.ip together with (jl.4p for sufficiently small e. Above that threshold, the 
droplet density of energy-minimizers converges to the uniform density jl. 

The key point that enables the analysis above is a kind of T-equivalence between the 
energy functional in (|l.lj) and its screened sharp interface analog (for general notions of 
T-equivalence or variational equivalence, see [2l[8]): 

E £ [u} = - I \Vu\dx + - [ I [u(x)-u £ )G(x-y)(u(y)-u e )dxdy. (1.9) 
2 h\ 2 h\ h\ 

Here, G is the screened potential as in (jl.6p . and u G A, where 

A:=BV(T 2 ;{-1,1}), (1.10) 

and we note that on the level of E £ the neutrality condition in (ll.3p has been removed. 
As we showed in [19], following the approach of [26], for £ £ given by (jl.ip in which u = u £ 
and u £ is defined in (ll.4p . we have 

minE e = rniuE 6 + 0(e a minE £ ), (1.11) 

for some a > 0. Therefore, in order to understand the leading order asymptotic expansion 
of the minimal energy minf £ in terms of | lne| _1 , it is sufficient to obtain such an expansion 
for mmE £ . This is precisely what we will do in the present paper. 
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In view of the discussion above, in this paper we concentrate our efforts on the analysis 
of the sharp interface energy E £ in (jl.9|) . An extension of our results to the original diffuse 
interface energy £ £ would lead to further technical complications that lie beyond the scope 
of the present paper and will be treated elsewhere. Here we wish to extract the next order 
non-trivial term in the T-expansion of the sharp interface energy E £ after (jl.5p . In contrast 
to [26], we will not use the Euler-Lagrange equation associated to (jl.9p . so our results about 
minimizers will also be valid for "almost minimizers" (cf. Theorem [2]). 

We recall that for e <C 1 the energy minimizers for E £ and 5 > 5 C consist of 0( \ lne|) 
nearly circular droplets of radius r ~ 3 1 / 3 e 1 / 3 | lne| -1 / 3 uniformly distributed throughout 
Tj [26l Theorem 2.2]. This is in contrast with the study of [121113] for a closely related 
energy, where the number of droplets remains bounded as e — > 0, and the authors extract 
a limiting interaction energy for a finite number of points. 

By T-convergence, we obtained in [191 Theorem 1] the convergence of the droplet density 
of almost minimizers (u £ ) of E £ : 

fjf( x ) : = ^" 2/3 | lne|- 1 / 3 (l + u £ (x)), (1.12) 

to the uniform density p, defined in (|1.8p . However, this result does not say anything about 
the microscopic placement of droplets in the limit e — > 0. In order to understand the 
asymptotic arrangement of droplets in an energy minimizer, our plan is to blow-up the 
coordinates by a factor of yj\ lne|, which is the inverse of the scale of the typical inter- 
droplet distance, and to extract the next order term in the T-expansion of the energy in 
terms of the limits as e — > of the blown-up configurations (which will consist of an infinite 
number of point charges in the plane with identical charge). 

We will show that the arrangement of the limit point configurations is governed by 
the Coulombic renormalized energy W, which was introduced in |34j . That energy W 
was already derived as a next order T-limit for the magnetic Ginzburg-Landau model of 
superconductivity [341135] . and also for two-dimensional Coulomb gases [37] . Our results 
here follow the same method of [35], and yield almost identical conclusions. 

The "Coulombic renormalized energy" is a way of computing a total Coulomb inter- 
action between an infinite number of point charges in the plane, neutralized by a uniform 
background charge (for more details see Section [2]). It is shown in [35] that its minimum is 
achieved. It is also shown there that the minimum among simple lattice patterns (of fixed 
volume) is uniquely achieved by the triangular lattice (for a closely related result, see [TU]). 
and it is conjectured that the triangular lattice is also a global minimizer. This triangular 
lattice is called "Abrikosov lattice" in the context of superconductivity and is observed in 
experiments in superconductors |41| . 

The next order limit of E £ that we shall derive below is in fact the average of the 
energy W over all limits of blown-up configurations (i.e. average with respect to the 
blow up center). Our result says that limits of blow-ups of (almost) minimizers should 
minimize this average of W. This permits one to distinguish between different patterns 
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at the microscopic scale and it leads, in view of the conjecture above, to expecting to see 
triangular lattices of droplets (in the limit e — > 0), around almost every blow-up center 
(possibly with defects). Note that the selection of triangular lattices was also considered in 
the context of the Ohta-Kawasaki energy by Chen and Oshita [10] , but there they were only 
obtained as minimizers among simple lattice configurations consisting of non-overlapping 
ideally circular droplets. 

It is somewhat expected that minimizers of the Ohta-Kawasaki energy in the macro- 
scopic setting are periodic patterns in all space dimensions (in fact in the original paper [28] 
only periodic patterns are considered as candidates for minimizers). This fact has never 
been proved rigorously, except in one dimension by Miiller [23] (see also [31LI42j ) . and at 
the moment seems very difficult. For higher-dimensional problems, some recent results in 
this direction were obtained in [2, 26, 38j establishing equidistribution of energy in various 
versions of the Ohta-Kawasaki model on macroscopically large domains. Several other 
results |12 l ll3 l[T5ll39j were also obtained to characterize the geometry of minimizers on 
smaller domains. The results we obtain here, in the regime of small volume fraction and in 
dimension two, provide more quantitative and qualitative information (since we are able to 
distinguish between the cost of various patterns, and have an idea of what the minimizers 
should be like) and a first setting where periodicity can be expected to be proved. 

The Ohta-Kawasaki setting differs from that of the magnetic Ginzburg-Landau model 
in the fact that the droplet "charges" (i.e., their volume) are all positive, in contrast with 
the vortex degrees in Ginzburg-Landau, which play an analogous role and can be both 
positive and negative integers. It also differs in the fact that the droplet volumes are not 
quantized, contrary to the degrees in the Ginzburg-Landau model. This creates difficulties 
and the major difference in the proofs. In particular we have to account for the possibility 
of many very small droplets, and we have to show that the isoperimetric terms in the energy 
suffice to force (almost) all the droplets to be round and of fixed volume. This has to be 
done at the same time as the lower bound for the other term in the energy, for example an 
adapted "ball construction" for non-quantized quantities has to be re-implemented, and 
the interplay between these two effects turns out to be delicate. 

Our paper is organized as follows. In Section [2] we formulate the problem and state 
our main results concerning the T-limit of the next order term in the energy (jl.9j) after 
the zeroth order energy derived in [19] is subtracted off. In Section [31 we derive a lower 
bound on this next order energy via an energy expansion as done in [19] however isolating 
lower order terms obtained via the process. We then proceed via a ball construction as 
in |20 P331I34] to obtain lower bounds on this energy in Section [4] and consequently obtain an 
energy density bounded from below with almost the same energy via energy displacement 
as in [35] in Section [5j In Section [6] we obtain explicit lower bounds on this density on 
bounded sets in the plane in terms of the renormalized energy for a finite number of points. 
We are then in the appropriate setting to apply the multiparameter ergodic theorem as 
in [35] to extend the lower bounds obtained to global bounds, which we present at the 
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end of Section Finally the corresponding upper bound (cf. Part (ii) of Theorem [I]) is 
presented in Section [7J 

Some notations. We use the notation (u £ ) £ A to denote sequences of functions u £ £ A 
as e = e n — > 0, where A is an admissible class. We also use the notation jjl € 
to denote a positive finite Radon measure d\x on the domain O. With a slight abuse of 
notation, we will often speak of /i as the "density" on Q and set dfi(x) = fi(x)dx whenever 
fi € L 1 (il). With some more abuse of notation, for a measurable set E we use \E\ to 
denote its Lebesgue measure, \dE\ to denote its perimeter (in the sense of De Giorgi), 
and n(E) to denote f E dfi. The symbols H 1 ^), BV(U), C k {n) and iJ-^fi) denote the 
usual Sobolev space, the space of functions of bounded variation, the space of fe-times 
continuously differentiable functions, and the dual of H 1 ^), respectively. The symbol 
o e (l) stands for the quantities that tend to zero as e — > with the rate of convergence 
depending only on £, 5 and k. 



2 Problem formulation and main results 

In the following, we fix the parameters K > 0, 5 > and £ > 0, and work with the energy 
E e in (jl.9p . which can be equivalently rewritten in terms of the connected components 
fif of the family of sets of finite perimeter f2 e := {u e = +1}, where (u £ ) £ A are almost 
minimizers of E £ , for sufficiently small e (cf. the discussion at the beginning of Sec. 2 
in pj5]). The sets Q £ can be decomposed into countable unions of connected disjoint sets, 
i.e., VL £ = Ui^!> whose boundaries dClf are rectifiable and can be decomposed (up to 
negligible sets) into countable unions of disjoint simple closed curves. Then the density \x e 
in (|1.12p can be rewritten as 

^( x ): =e -2/3| lne |-i/3^ x ^ (;c)5 (21) 

i 

where xn? are the characteristic functions of fif . Motivated by the scaling analysis in the 
discussion preceding equation (|1.12p . we define the rescaled areas and perimeters of the 
droplets: 

A\ := e~ 2/3 \ lne| 2 / 3 |Of |, Pf := e~ 1/3 \ lne\ 1/3 \dnf\. (2.2) 

Using these definitions, we obtain (see }19ll26j) the following equivalent definition of the 
energy of the family (u £ ): 

E £ [u £ ] = e 4 / 3 | lne| 2 / 3 + E £ [u £ ]\ , (2.3) 
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where 



E £ [u £ ] := E (*? ~ % A i) + 2 // G ( x ~ vW(x)dtf(y). (2.4) 

Also note the relation 

M £ (Tl) = n ^E^- ( 2 - 5 ) 

i 

As was shown in |19|J26]. in the limit e — > the minimizers of i? £ are non-trivial if and 
only if 5 > 5 C , and we have asymptotically 

mmE £ ~ -4r(25 - 5 c )e 4/3 \ lne\ 2/3 £ 2 as e -> 0. (2.6) 

Furthermore, if /U £ is as in (|2.ip and we let t> £ be the unique solution of 

-Av £ + k 2 v £ = ^ in W 2 ' p (T 2 e ), (2.7) 
for any p < oo, then we have 

« e ^:=-^-*c) m ^(T, 2 ). (2.8) 

To extract the next order terms in the T-expansion of E e we, therefore, subtract this 
contribution from E e to define a new rescaled energy F £ (per unit area): 

F £ [u] := e" 4/3 1 In e\ 1/3 r 2 E £ [u] - I In el (25 - 4) + — 3— (5 - 5 C ) (hi I In el + In 9) . (2.9) 

Note that we also added the third term into the bracket in the right-hand side of (|2,9p 
to subtract the next-to-leading order contribution of the droplet self-energy, and we have 
scaled F £ in a way that allows to extract a non-trivial 0(1) contribution to the minimal 
energy (see details in Section [3]). The main result of this paper in fact is to establish 
T-convergence of F £ to the renormalized energy W which we now define. 

In [35], the renormalized energy W was introduced and defined in terms of the su- 
perconducting current j, which is particularly convenient for the studies of the magnetic 
Ginzburg-Landau model of superconductivity. Here, instead, we give an equivalent def- 
inition, which is expressed in terms of the limiting electrostatic potential of the charged 
droplets, after blow-up, which is the limit of some proper rescaling of v £ (see below). How- 
ever, this limiting electrostatic potential will only be known up to additive constants, due 
to the fact that we will take limits over larger and larger tori. This issue can be dealt with 
in a natural way by considering equivalence classes of potentials, whereby two potentials 
differing by a constant are not distinguished: 

[cp] ■= {cp + c | cG E}. (2.10) 
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This definition turns the homogeneous spaces W rl -P(M d ) into Banach spaces of equivalence 
classes of functions in W^(M. d ) defined in (|2.1(jp (see, e.g., [25]). Here we similarly define 
the local analog of the homogeneous Sobolev spaces as 



):={H \<peWk P c (R 2 )}, (2.11) 



with the notion of convergence to be that of the L^ oc convergence of gradients. In the 
following, we will omit the brackets in [•] to simplify the notation and will write ip G 
W^(R 2 ) to imply that (p is any member of the equivalence class in (|2,10D , 
We define the admissible class of the renormalized energy as follows : 



Definition 2.1. For given m > and p G (1,2), we say that ip belongs to the admissible 



class Am, if '</? G Whf{M?) and ip solves distributionally 



-Aip = 2vr^5 a - m, (2.12) 

aeA 

where AcM 2 is a discrete set and 

lim -|r / 5 a (x)dx = m. (2.13) 

R^oc R* j Br{0) ^ 

Remark 2.2. Observe that if (p G A m , then for every x G B R (0) we have 

<p(x) = ln \ x - + m(x), (2-14) 

where A R := An B R (0) is a finite set of distinct points and ip R G C°°(IR 2 ) is analytic in 
B R (0). In particular, the definition of A m is independent of p. 

We next define the renormalized energy. 
Definition 2.3. For a given ip G A m , the renormalized energy W of tp is defined as 

m>0 

W(tp) :=limsuplim— ^ ( / i| Vip\ 2 X Rdx + vr In 77 V xr(p) J > i 2 - 1 ^ 

v^o \K R \ \Jm?\U aeA B n (a) 2 tth ) 

where K R = [— R, R] 2 , \R is a smooth cutoff function with the properties that < \R < 1) 
in K R \ (dK R U K R _i), Xr(%) = 1 f or al1 x € Kr-\, Xr(%) = for all x G R 2 \K R , and 
|Vxh| < C for some C > independent of R. 



Various properties of W are established in |35| . we refer the reader to that paper. The 
most relevant to us here are 
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1. min_4 m W is achieved for each m > 0. 

2. If ip G «4. m and <^'(x) := </?(-^=), then y>' E and 

W(p) = m \W{ip') - ^ log m 

hence 

min W = m [ min W — - log m 

A m \Ai 4 

3. is minimized over potentials in A\ generated by charge configurations A consisting 
of simple lattices by the potential of a triangular lattice, i.e. |35^ Theorem 2 and 
Remark 1.5], 

min W(ip) = W(ip A ) = -- ln(V2^b \rj(r)\ 2 ) ~ -0.2011, 

<p£A\ 2 
A simple lattice 

where r = a + ib, rj(r) = g 1//24 r[n>i(l ~~ 9™) 1S the Dedekind eta function, q = e 2nlT , 

a and b are real numbers such that = ^(1, 0)Z ® (a, b)7^j is the dual lattice 

to a triangular lattice A A whose unit cell has area 2ir, and y? A solves (|2.12j) with 
A = A A . 

In particular, from property 2 above it is easy to see that the role of m in the definition of 
W is inconsequential. 

We are now ready to state our main result. Let £ £ := | lnel 1 / 2 ^. For a given u £ £ A, 
we then introduce the potential (recall that ip £ is a representative in the equivalence class 
defined in (I2TT0]) ) 

ip £ (x) := 2 • 3~ 2/3 | lne| v £ (x\ lne\~ 1/2 ), (2.17) 

where v £ is a periodic extension of v £ from T| e to the whole of K 2 . We also define V to be 
the family of translation-invariant probability measures on Wi (M 2 ) concentrated on ,4m 
with m = 3" 2 / 3 (<5- <5 C ). 

Theorem 1. (r-convergence of F £ ) Fix k > 0, 5 > 5 C , p € (1, 2) and £ > 0, and Zei F e 

6e defined by S2. 9\) . Then, as e — > we have 



F £ 4 F°[P] := 3 4 / 3 y W( V )dP(^) + 32/3( g * c) , (2.18) 
where P £ V . More "precisely: 
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i) (Lower Bound) Let (u £ ) £ A be such that 

limsup F £ [u £ ] < +00, (2.19) 

and let P 6 be the probability measure on W,jf (M 2 ) which is the pushforward of the 
normalized uniform measure on T 2 e by the map x i-> ip £ (x+-), where ip £ is as in (|2.17p . 
Then, upon extraction of a subsequence, (P £ ) converges weakly to some P G V , in 
the sense of measures on W^(M. 2 ) and 

liminf F e [u £ ] > F 



(2.20) 



ii) (Upper Bound) Conversely, for any probability measure P G V , letting Q be its push- 
forward under —A, there exists {u £ ) G A such that letting Q £ be the pushforward 
of the normalized Lebesgue measure on T| e by x ^ —A(p £ (x + •), where (p £ is as in 
(|2.17p . roe Ziaue Q £ — v Q, in i/te sense 0/ measures on W^' P (M. 2 ), and 

limsupF e [n £ ] < F°[P], (2.21) 

as e — > 0. 

We will prove that the minimum of F° is achieved. Moreover, it is achieved for any P G V 
which is concentrated on minimizers of A m with m = 3~ 2 / 3 (5 — <5 C ). 

Remark 2.4. The phrasing of the theorem does not exactly fit the framework of T -convergence, 
since the lower bound result and the upper bound result are not expressed with the same 
notion of convergence. However, since weak convergence of P e to P implies weak con- 
vergence of Q e to Q, the theorem implies a result of Y -convergence where the sense of 
convergence from P £ to P is taken to be the weak convergence of their push-forwards Q £ to 
the corresponding Q. 

The next theorem expresses the consequence of Theorem [T] for almost minimizers: 

Theorem 2. Let m = 3 _2 / 3 (<5 — <5 C ) and let {u £ ) £ A be a family of almost minimizers of 
F°, i.e., let 

lim F £ \u £ ] = minF . 
e^o v 

Then, if P is the limit measure from Theorem^ P -almost every ip minimizes W over A m . 
In addition 

rninF = 3 4 / 3 minl^ + ~ . (2.22) 

V Am 8 ^ ' 
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Note that the formula in (|2.22p is not totally obvious, since the probability measure con- 
centrated on a single minimizer ip S A m of W does not belong to V . 

The result in Theorem [2] allows us to establish the expansion of the minimal value of 
the original energy E £ by combining it with (|2.9|) and (jl.lip . 

Theorem 3. (Asymptotic expansion of min£ e ) Let V = ^(1 — u 2 ) 2 , k = | and 
m = 3~ 2//3 (<5 — 5 C ). Fix 5 > 5 C and I > 0, and let £ £ be defined by with u = u £ from 

(TTT^p. Then, as e —¥ we have 

r 2 mm£ £ = -^-(25-S c )e i/3 \lne\ 2/3 - ~5 c )e i/3 \ In el" 1 / 3 (In | lnel + In 9) 

+e 4 / 3 \lne\- 1 / 3 ( 3 4 / 3 min W + ~ ) + o(e 4 / 3 | ln^" 1 / 3 ). 

\ Am 8 / 

(2.23) 

As mentioned above, the T-limit in Theorem [1] cannot be expressed in terms of a single 
limiting function (p, but rather it effectively averages W over all the blown- up limits of <p £ , 
with respect to all the possible blow-up centers. Consequently, for almost minimizers of the 
energy, we cannot guarantee that each blown-up potential ip £ converges to a minimizer of 
W, but only that this is true after blow-up except around points that belong to a set with 
asymptotically vanishing volume fraction. Indeed, one could easily imagine a configuration 
with some small regions where the configuration does not ressemble any minimizer of W, 
and this would not contradict the fact of being an almost minimizer since these regions 
would contribute only a negligible fraction to the energy. Near all the good blow-up centers, 
we will know some more about the droplets: it will be shown in Theorem [J] that they are 
asymptotically round and of optimal radii. 

We finish this section with a short sketch of the proof. Most of the proof consists in 
proving the lower bound, i.e. Part (i) of Theorem[TJ The first step, accomplished in Section 
[3] is, following the ideas of [26], to extract from F £ some positive terms involving the sizes 
and shapes of the droplets and which are minimized by round droplets of fixed appropriate 
radius. These positive terms, gathered in what will be called M e , can be put aside and 
will serve to control the discrepancy between the droplets and the ideal round droplets 
of optimal sizes. We then consider what remains when this M e is subtracted off from F £ 
and express it in blown-up coordinates x' = x*J\ lne|. It is then an energy functional, 
expressed in terms of some rescaling of (p £ which has no sign and which ressembles that 
studied in [35]. Thus we apply to it the strategy of [35]. The main point is to show that, 
even though the energy density is not bounded below, it can be transformed into one that 
is by absorbing the negative terms into positive terms in the energy in the sense of energy 
displacement [35], while making only a small error. In order to prove that this is possible, 
we first need to establish sharp lower bounds for the energy carried by the droplets (with 
an error o(l) per droplet). These lower bounds contain possible errors which will later be 
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controlled via the M e term. This is done in Section[3]via a ball construction as in [20|83p34j. 
In Section we use these lower bounds to perform the energy displacement as in [35]. Once 
the energy density has been replaced this way by an essentially equivalent energy density 
which is bounded below, we can apply the abstract scheme of [35] that serves to obtain 
lower bounds for "two-scale energies" which T-converge at the microscopic scale, via the 
multiparameter ergodic theorem. This is achieved is Section [6j Prior to this we obtain 
explicit lower bounds at the microscopic scale in terms of the renormalized energy for a 
finite number of points. It is then these lower bounds that get integrated out, or averaged 
out at the macroscopic scale to provide a global lower bound. 

Finally, there remains to obtain the corresponding upper bound. This is done via an 
explicit construction of a periodic test-configuration, following again the method of [35] . 

3 Derivation of the leading order energy 

In preparation for the proof of Theorem [H we define 

p e :=3 1 / 3 e 1 / 3 |ln e | 1 /6 and f e :=(j^X\ (3.1) 

VI 111 Pe\J 

Recall that to leading order the droplets are expected to be circular with radius S 1 / 3 ^ 1 / 3 1 In e| 
Thus p £ is the expected radius, once we have blown up coordinates by the factor of yf\ lne|, 
which will be done below. Also, we know that the expected normalized area M is 3 2 / 3 
but this is only true up to lower order terms which were negligible in j!9j ; as we show below, 
a more precise estimate is Ai ~ nf £ , so f e above can be viewed as a "corrected" normal- 
ized droplet radius. Since our estimates must be accurate up to o e (l) per droplet and the 
self-energy of a droplet is of order Af lnp e , we can no longer ignore these corrections. 

The goal of the next subsection is to obtain an explicit lower bound for F e defined by 
(|2.9p in terms of the droplet areas and perimeters, which will then be studied in Sections 
|4]and onward. We follow the analysis of |19j . but isolate higher order terms. 



3.1 Energy extraction 

We begin with the original energy E £ (cf. (|2.4j) ) while adding and subtracting the truncated 
self interaction: first we define, for 7 € (0, 1), truncated droplet volumes by 

7. . = / A ? if AI < 3 2 / 3 vr 7 -\ 

*' {(S^ir^lAfl) 1 / 2 if Af > 32/3,^-1, { - } 

as in [19] . The motivation for this truncation will become clear in the proof of Proposi- 
tion when we obtain lower bounds on the energy on annuli. In [TU] the self-interaction 

energy of each droplet extracted from E e was ^-1 ^ 1 , yielding in the end the leading order 
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energy E° [fj] in (jl.5p . A more precise calculation of the self-interaction energy corrects the 
coefficient of \Af\ 2 by an 0(ln | lne|/| lne|) term, yielding the following corrected leading 
order energy for E £ : 

E °elt4 ■= CI + fr " / + 2 [[ G{x- y)diM(z)Mv)- (3-3) 

The energy in (j3.3f) is explicitly minimized by df/,(x) = p, E dx (again a correction to the 
previously known p, from (jl.8p ) where 

1 3k 2 \ , - 3k 2 .„ 



and 



mi "^=^l 25 (l) 1/3 -^(l) 2/ T < 35 > 



Observing that f e — >■ 3 1//3 we immediately check that 

fi £ — > p as e — > 0, (3-6) 

and in addition that (|3.5p converges to the second expression in (|1.8j) . To obtain the next 
order term, we Taylor-expand the obtained formulas upon substituting the definition of f £ . 
After some algebra, we obtain 

„ 9 <5c /„r f\ 1/1= -f s In I In el + ln9 _ / (In I lne|) 2 \ .„ „. 

r" m mE ? = ^ (M - 4) - j^C* - <J ^ + O (L-LJL j . (3.7) 

Recalling once again the definition of F e from (j2.9j) . we then find 

F £ [u £ ] = |lne| (e~ 4 / 3 | lne| ~ 2 / 3 r 2 i? £ [u £ ] - r 2 min^ E °) + O ( ^r^p 

and in view of the definition of -E £ from (|2.3p . we thus may write 

^ 2 / (In |lne 

2k 2 V ' " V line 

Thus obtaining a lower bound for the first term in the right-hand side of (j3.8j) implies, up 
to o e (l), a lower bound for F 5 . This is how we proceed to prove Lemma 13. II below. 
With this in mind, we begin by setting 

K 1 (r 3k 2 \ . , 

v ' = * + u^' v =2^ [ s -xr < (3 - 9) 



F £ [« £ ] = | lne|r 2 ( £ £ [u £ ] + -^2- - min£ £ ° j + O ( V 1— ^ ) . (3i 



where v £ is the solution to (|2.7p with right side equal to fi £ in (|3.4j) . 
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3.2 Blowup of coordinates 

We now rescale the domain Tj by making the change of variables 



x = 



VI me l> 

h' £ (x') = h £ (x), (3.10) 



Observe that 



i £ = £y/\kie\. 

tp £ (x') = 2 • 3- 2 / 3 h' £ (x') Vx' G T|, (3.11) 



where (p £ is defined by (|2.17j) . It turns out to be more convenient to work with h' £ and 
rescale only at the end back to (p £ . 

3.3 Main result 

We are now ready to state the main result of this section, which provides an explicit lower 
bound on F £ . The strategy, in particular for dealing with droplets that are too small or 
too large is the same as [19] . except that we need to go to higher order terms. 



Proposition 3.1. There exist universal constants 7 G (0, g), c\ > 0, C2 > 0, C3 > and 
£q > such that if 5 > 5 C and (u £ ) G A with £l £ := {u £ > 0}, then for all e < £q 

eF £ [u £ ] >M £ + -^- [ (\Vh'f + j^lKl 2 ) dx> - 4j £ |if | 2 + 0.(1), 
me 7 T 2 V hie / nrf ^ 

(3.12) 



Af>3 2 / 3 7T7 



where M £ > is defined by 
M -=E(^-V / 4^4f)+ c i E A * 

i Af>3 2 / 3 7T7- 1 

+ c 2 £ (Af - nflf + c 3 £ Af. (3.13) 

32/3 vr7 <yl|<32/3 7r7 -l Af<3 2 / 3 7r 7 

Remark 3.2. Defining j3 := 3 2//3 7T7, fry isoperimetric inequality applied to each connected 
component of £l £ separately every term in the first sum in the definition of M £ in (|3.13p 
is non-negative. In particular, M £ measures the discrepancy between the droplets Of with 
Af > (3 and disks of radius f £ . 
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The proposition will be proved below, but before let us examine some of its fur- 
ther consequences. The result of the proposition implies that our a priori assumption 
lim sup e _ 5>0 F e [u £ ] < +00 translates into 



\lne\ J T 2 V |ln £ |' J Trrf 1 



for some C > independent of e <C 1, which, in view of (|3.ip is also 

A major goal of the next sections is to obtain the following estimate 

W7l (L H 2 ^"^'-^""*! E l*l 2 ) * -Cln*(M. + 2), (3.15) 

for some C > independent of e <C 1, so that the a priori bound (|3,14p in fact implies 
that M e is uniformly bounded independently of e for small e. This will be used crucially 
in Section HOI 

We note that h' e (x') satisfies the equation 



where we define in T 2 £ 



and 



-Ah' e + -^-h' e = ^-^ mW 2 *(T%) (3.16) 



£&):=Y,^W), (3-17) 



SKx') := (3.18) 

which will be used in what follows. Notice that each Sf(x') approximates the Dirac delta 
concentrated on some point in the support of ft' i£ and, hence, fx' E (x')dx' approximates the 
measure associated with the collection of point charges with magnitude A\. In particular, 
the measure d\J e evaluated over the whole torus equals the total charge: ^4(T 2 £ ) = A\. 
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3.4 Proof of Proposition 13.11 

- Step 1: We are first going to show that for universally small e > and all 7 G (0, g) we 
have 

£ 2 F £ [u e ] > Ti + T 2 + T 3 + T 4 + T 5 + o e (l), (3.19) 



where 



T i = E " V 7 ^! ) > (3-20) 



//2 



T 2 = l^ E (^"^1) 2 > (3-21) 

3 2 / 3 7T7</l?<3 2 / 3 7r7- 1 

T 3 = ^^ E ^f(^--1) 2 > (3-22) 

Af<3 2 / 3 7T7 

T 4 = 2 (6-V 1 -!)^, (3.23) 

Af>3 2 / 3 7r 7 - 1 

T 2 
5 Unel 



I (\Vh £ \ 2 + K 2 \h e f)dx--^ V|if| 2 . (3.24) 

To bound -F^it 5 ] from below, we start from (|3.8p . In particular, in view of (|2.7p we may 
rewrite (|2.4p as 

^K] = _L-Wi?_ Sa £ ) + 2 / (|v^| 2 + *V| 2 ) rfx 

i 

+ 2 / (^Wl^x--^ £|if| 2 . (3.26) 

We start by focusing on (|3.25p . First, in the case A\ > 3 2 / 3 7T7 _1 we have \Af\ 2 = 

3 2//3 7T7 _1 J 4? and hence, recalling that f e = 3 1 / 3 + o e (l), where o e (l) depends only on 
e, we have for e universally small and 7 < i: 



7rr| 7rr| 



>^(I + i( 7 -i-6)). (3.27) 
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We conclude that for A\ > 3 2 / 3 7T7 1 , we have 

1£|2 



(3.28) 



On the other hand, when < 3 2//3 7T7 1 we have = A? and we proceed as follows. 
Let us begin by defining, similarly to [19], the function 



yx 7rr| 

for x € (0, +oo) and observe that / is convex and attains its minimum of J- at x = 7rf 2 , 
with 

™ = ^>°- 

By a second order Taylor expansion of / around 7rr 2 , using the fact that /" is decreasing 
on (0, +oo), we then have for all x < xq 

v / 4^+ ^ = xf{x) > x I A + ^ (x - vrf 2 ) 2 ] . (3.29) 

nr e \ r £ AXq J 

We, hence, conclude that when 3 2 / 3 7T7 < A\ < 3 2 / 3 7T7 _1 , we have 

v 71 ^ + %- 2 i Al -{v £ - 2 i) Al + ^^ AKAt ~ 7rr1)2 ' (3 - 30) 

and when Af < 3 2/,3 7T7, we have 

Combining (|3.28[> . (|3.30j) and (|3.31j) . summing over all i, and distinguishing the different 
cases, we can now bound (|3.25p from below as follows: 



7 /2 

3 2 / 3 7r 7 <ylf<3 2 / 3 7r7- 1 
7 -5/2 

+ 4yr 2 . 3 2/3 Z> ^i(^i-^e) 

Af<3 2 / 3 7T7 

+ E (6-V 1 -!)^!- (3-32) 

Af>3 2 / 3 7T7- 1 
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We now focus on the term in (|3.26p . Using (|3.9p . we can write the integral in (|3.26j) as: 

2 I (Vv £ \ 2 + K 2 \v £ \ 2 )dx 

= —?—[ (\Vh £ \ 2 + K 2 h 2 ) dx [ h £ dx + 2K 2 \v e \ 2 e. (3.33) 

|lne| 2 Jtj v £J |hie| h\ 

Integrating (|2.7p over T 2 and recalling the definition of h £ in (|3.9p . as well as (|2.5p . leads 

y 2 /i £( ix = ^^^^f-4 K > e | 2 ^. (3.34) 



to 

4k 2 v £ f , , 4v £ 



| hie 

Combining (|3.33p and (|3.34[) . we then find 

2/ (|Vv £ | 2 + K> e | 2 ) dx = — [ (\Vh e \ 2 + K 2 h 2 ) dx 
Jrj I hie I Jjj 

Also, by direct computation using (13. 5|) and (13. 9p we have 

2K 2 \v £ \ 2 £ 2 = ^-mmE° £ . (3.36) 



Therefore, combining this with (|3.8|) . (|3.32|) and (|3.35|) . after passing to the rescaled coor- 
dinates and performing the cancellations we find that 

£ 2 F £ [u £ ] >T 1 + T 2 + T 3 + T 4 + -^- [ ( \Vh' £ (x')\ 2 + J^-\h' £ (x')\ 2 ) dx' 

|lne| Jf2 e \ |lne| J 

-^a£W+*(i). ( 3 - 37 ) 

c i 

which is nothing but (|3.19p . 

- Step 2: We proceed to absorbing the contributions of the small droplets in (|3.24p by 
( 13.211) . To that effect, we observe that, for the function 

^ ■= -^wA x - - - r r1 ~ ( 27rr1 + ) 4 ' (3 ' 38) 



there exists a universal 7 G (0, |) such that <& £ (%) > x whenever < x < 3 2 ^ 3 ir-y and e is 
universally small. Using this observation, we may absorb all the terms with Af < 3 2 / 3 7T7 
appearing in the second term in f|3.24|) into (|3.22[) by suitably reducing the coefficient in 
front of the latter. This proves the result. □ 
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4 Ball construction 



The goal of this section is to show (|3.15|) using the abstract framework of Theorem 3 in [35] . 
The difficulty in doing this, as in the case of the Ginzburg-Landau model treated in [35], 
is that the energy density e' £ — ^| hip e \ Xm e >/3 l^f l 2< ^f ^ s no * positive (or bounded below 
independently of (u £ )). The next two subsections are meant to go around this difficulty 
by showing that this energy density can be modified, by displacing a part of the energy 
from the regions where the energy density is positive into regions where the energy density 
is negative in order to bound the modified energy density from below while making only 
a small enough error. This is achieved by obtaining sharp lower bounds on the energy of 
the droplets. Since their volumes and shapes are a priori unknown, the terms in M E are 
used to control in a quantitative way the deviations from the droplets being balls of fixed 
volume. 

In this section we perform a ball construction which follows the procedure of [35]. The 
goal is to cover the droplets {1^ e } whose volumes are bounded from below by a given 
(3 > with a finite collection of disjoint closed balls whose radii are smaller than 1, on 
which we have a good lower bound for the energy in the left-hand side of (I3.15|) . This 
is possible for sufficiently small e in view of the fact that £ e — > oo and that the leading 
order asymptotic behavior of the energy from f|2.6[) yields control on the perimeter and, 
therefore, the essential diameter of each of £l' ie - The precise statements are given below. 
We will also need the following basic result, which holds for sufficiently small e ensuring 
that the droplets are smaller than the sidelength of the torus (see the discussion at the 
beginning of Sec. 2 in |19|). 

Lemma 4.1. There exists Eq > depending only on i, k, 5 and sup £>0 F e [u £ ] such that 
for all e < Eq we have 

essdiam(^ jE )<c|c>ny, (4.1) 

for some universal c > 0. 

^From now on and for the rest of the paper we fix 7 to be the constant given in 
Proposition 13.11 and, as in the previous section, we define (3 = 3 2 / 3 irj. We also introduce 
the following notation which will be used repeatedly below. To index the droplets, we will 
use the following definitions: 

Ip := {i G N : A\ > /?}, I E := {i G N : D {T%\E)\ = 0}, I PjE := I p D I E , (4.2) 

where E C T| e . For a collection of balls B, the number r(B) (also called the total radius 
of the collection) denotes the sum of the radii of the balls in B. For simplicity, we will say 
that a ball B covers Vt^ e , if i G 1b- 

The principle of the ball construction introduced by Jerrard [20] and Sandier [33| and 
adapted to the present situation is to start from an initial set, here (Jig/ v e ^ or a gi ven 
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U C and cover it by a union of finitely many closed balls with sufficiently small radii. 
This collection can then be transformed into a collection of disjoint closed balls by the 
procedure, whereby every pair of intersecting balls is replaced by a larger ball whose radius 
equals the sum of the radii of the smaller balls and which contains the smaller balls. This 
process is repeated until all the balls are disjoint. The obtained collection will be denoted 
Bo, its total radius is r(Bo)- Then each ball is dilated by the same factor with respect to 
its corresponding center. As the dilation factor increases, some balls may touch. If that 
happens, the above procedure of ball merging is applied again to obtain a new collection of 
disjoint balls of the same total radius. The construction can be stopped when any desired 
total radius r is reached, provided that r is universally small compared to t e . This yields 
a collection B r covering the initial set and containing a logarithmic energy |2Uy33j. 

We now give the statement of our result concerning the ball construction and the 
associated lower bounds. Throughout the rest of the paper we use the notation / + := 
max(/, 0) and /~ := — min(/, 0). 

Proposition 4.2. Let U C T^ e be an open set such that Ip y ^ 0, and assume that (12.19|) 
holds. 

- There exists Eq > 0, ro £ (0, 1) and C > depending only on £, k, 6 and sup £>0 -F £ [u £ ] 
such that for all e < Eq there exists a collection of finitely many disjoint closed balls 
Bo whose union covers \J ie j g v £ and such that 

r(Bo) < ce 1/3 \ lne\ 1/6 ^ Pf < r , (4.3) 

for some universal c > 0. Furthermore, for every r £ [r(Bo),ro] there is a family of 
disjoint closed balls B r of total radius r covering Bo- 

- For every B £ B r such that B <ZU we have 

[ (\Vh' £ \ 2 dx' + -^-\h' £ \ 2 )dx'> — (\n^—-cr) + V |If| 2 , 
JbV d 4|lne|' d ; " 2tt V r(B ) J ^ 1 ll 

for some c > depending only on k and 5. 

- If B G B r , for any non-negative Lip schitz function x with support in U, we have 

> -c\\Vx\\oo E l^l 2 ' 

where Xi '■= JuX^id x '; with 5f(x') defined in (|3.18p . for some c > depending only 
on k and 6, and a universal C > 0. 
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Remark 4.3. The explanation for the factor of | in front of p^l^el 2 is that we must 
'save ' a fraction of this term for the mass displacement argument in Section [5| and in the 
convergence result in Section® 

Proof of the first item. Choose an arbitrary ro € (0, 1). As in [19] . from the basic lower 
bound on E £ (see [El Equations (2.12) and (2.15)]): 

*m >- ]t\ e i - j^Li e * + (e 2 • <«> 

where and are defined in (|2.2j) . we obtain with the help of (|2.19p that 

lim sup — *— V AS < C, lim sup - — r V Pf < C, (4.5) 
ho Ins V ho |lne|4^ 

2 



for some C > depending only on £, n, 5 and sup £>0 F £ [?x £ ]. 

As is well-known, the essential diameter of a connected component of a set of finite 
perimeter on a torus can be bounded by its perimeter, provided that the latter is universally 
small compared to the size of the torus (see, e.g., [1]). Therefore, in view of the the definition 
of Pf in (12, 2h and the second of (|4.5p , for sufficiently small e it is possible to cover each £ 
with i E Ip^j by a closed ball -Bj, so that the collection Bq consisting of all Bi's (possibly 
intersecting) has total radius 

roCBo) < Ce^llnel 1 / 6 Pf, (4.6) 

for some universal C > 0. Furthermore, by the first inequality in (14. 5p and the fact that 
Af > {3 for all z G J« j/ the collection Bo consists of only finitely many balls. Therefore, we 
can apply the construction a la Jerrard and Sandier outlined at the beginning of this section 
to obtain the desired family of balls Bq and B r , with r(Bo) = t{Bq). The estimate on the 
radii follows by combining the second of (j4.5|) and (|4.6p and the fact that t e — > oo with the 
rate depending only on £, for sufficiently small e depending on £, At, 6, sup £>0 F e [u £ ] and 
r - □ 

Proof of the second item. Let B C U be a ball in the collection B r . Denote the radius of 
B by re and set 



K 2 



X F := —— I h'dx' 



£ 

B 



, l ne l 

Integrating (|3.16p over B and applying the divergence theorem, we have 



\ dH 1 (x')=m B ,e-X E , (4.7) 

J dB <JV 
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where 

m Bi£ := f - tfW = E ^ + E 9 ^ ~ ^l B l> 

Jb iei B i0 B 

for some 9{ £ [0, 1) representing the volume fraction in B of those droplets that are not 
covered completely by B, and v is the inward normal to dB. Using the Cauchy-Schwarz 
inequality, we then deduce from fj4.7|) that 

|V/4| 2 dU\x') > -L_ (m B>e - X £ ) 2 > (4.8) 

SB 27rrB 27rr B 

By another application of the Cauchy-Schwarz inequality, we may write 

f \ h 'f dx '> X {^ . (4.9) 
4| lne| J B 6 4irr B k 2 

We now add (|4.8p and (|4.9p and optimize the right-hand side over X e . We obtain 



i^r^ > ( i - (4.io) 

4|lne| V 2vrr B V I lne l7 

for C = fc . Recalling that r B < r < < 1, we can choose e sufficiently small depending 
only on k so that the term in parentheses above is positive. 

Inserting the definition of m B ^ £ into (|4,10p and discarding some positive terms yields 

2 



/ \Vti e \ 2 dH 1 {x') + -^— [ \h' £ \ 2 dx' 
JdB 4|lne| V 



^(E^ + E^f-^i) ! (-^) 

i£l B i0B 

>-^( E E ^) 2 (' - 2,r|B| (E " w I ■ < 411 » 

We now use the fact that by construction B covers at least one with Af > {3. This 
leads us to 
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for some c > depending only on k and 6, where in the last line we used that A\ > Af. 
Hence there exists ro £ (0, 1) depending only on k, and 5 such that the right-hand side of 
(|4.12p is positive. 

Finally, let us define T(x,r) := f B{x ^ \Vh' e \ 2 dx' + \h' £ \ 2 dx' , where B(x,r) 

is the ball centered at x of radius r. The relation (|4.12p then reads for B(x,r) = B £ B r 
and a.e. r £ (r(£>o),r ]: 

with c as before. Then using [341 Proposition 4.1], for every B £ $(s) := Z3 r with r = 
e s r(^o) (using the notation of |34| Theorem 4.2]) we have 

/.J^' + ^/j^^I E S E ATP— CW»* 

B'C-B 

= f E ^ E l^l 2 (l-«*r(Bb))(ft 

B'cB 

^Ei^ 2 ('"^)-). 

where we observed that the double summation appearing in the first and second lines is 
simply the summation over In b- Once again, in view of the fact that tb < 1 and that 
both terms in the integrand of the left-hand side of ()4.14|) are non-negative, this completes 
the proof of the second item. □ 

Proof of the third item. This follows [35] , Let x be a non-negative Lipschitz function with 
support in U. By the "layer-cake" theorem [21], for any B £ B r we have 

L* H 2 + ^' 2 ) « - CL ( |Vfc - 12 + ^k Kl2 ) (4 - l5) 

where E t := {x > t}. If i £ Ip.B, then by construction for any s £ [r(Bo),r] there exists a 
unique closed ball Bi jS £ B s containing Q! ie . Therefore, for t > we can define 

s(i,t) := sup{s £ [r(B ),r] : C £*} , 

with the convention that s(i,i) = r(Bo) if the set is empty. We also let B\ := B is r it \ 
whenever s(i,t) > r(Bo). Note that for each i £ Ip b we have that t *— > s(i,t) is a non- 
increasing function. In particular, we can define ti > to be the supremum of the set of 
i's at which s(i,t) = r (or zero, if this set is empty). 
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If t > ti and s(i,t) > r(£>o), then for any x £ Q' ie and any y G Bj\ E t (which is not 
empty) we have 

X (x) - t < X (x) - X (y) < 2«(i,t)||Vx||oo- (4.16) 
Averaging over all x G f2' i£ , we hence deduce 

Xi-t<2s{i,t)\\Vx\W (4.17) 

Now, for any t > the collection {l?*}^^ B t , where Ip,B,t '■= {i G I/3,B ■ s(i,t) > r(Bo)} 
is disjoint. Indeed if i,j G Ip,B,t and s(i,i) > s(j,t) then, since £> s (i,t) is disjoint, the balls 
Bi, s (i,t) an d are either equal or disjoint. If they are disjoint we note that s(i, t) > 

s(j,t) implies that £,>(j,t) C B jjS ^, and, therefore, B* = B j>s ^^ and J3| = -B M (i,t) 
are disjoint. If they are equal and s(i,t) > s(j,t), then Bj^ s ^^ C E t , contradicting the 
definition of s(j,t). So s(j,t) = s(i,t) and then B l - = B\. 

Now assume that B' G {S*}j g / j9Bt and let s be the common value of s(i,t) for z's in 
Ifi,B'- Then, the previous item of the proposition yields 



B' 



vh 'f + j^\ h 'A dx '>]_U n " cs ) y \a 

A\ In 6-1 1 61 I — 9-tt \ r(R n \ I ^ 1 

r 8,S',t 



4|lne|' £l J ~ 2ir \ r(B ) 



c |2 



Summing over B' G {fi|}j g j„ B1 , we deduce 

JbheA £ ' 4|lne| l£i ; " 2tt ^ 1 d V r(B„) 1 7 



+ 



r E i^fi 2 (in%4--(M)) , (4.18) 



where in the last inequality we took into consideration that all the terms corresponding 
to i G I/3 t B \ Ip,B,t give no contribution to the sum in the right-hand side. Integrating the 
above expression over t and using the fact that ro(Bo) < s(i,t) < r yields 



E t DB 



Vh >f + J^_\ h 'A dx ' dt> _L y TVln 414- cr) dt 

d 4 lne 1 £l J ~ 2ir ^ 1 d Jo V K^o ' 



We now concentrate on the last term in (|4.19p . Using the estimate in (|4.17p and the 
definition of ij, we can bound the integral in this term as follows 

^ ta fM*>r ta (_^_) (B >_ ||v x u, ,4,o) 
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for some universal C > 0, which is obtained by an explicit computation and the fact that 
r < ro < 1. Finally, combining (|4.'20p with (|4.19|) . the statement follows from (|4.15|) . □ 

Remark 4.4. Inspecting the proof, we note that the statements of the proposition are still 
true with the left-hand sides replaced by /mg xl^ 7 ^) 2 dx' + jrg^r j B xW £ \ 2 dx' (with % = 1 
or x Lipschitz, respectively). 



5 Energy displacement 



In this section, we follow the idea of [35] of localizing the ball construction and combine 
it with a "energy displacement" which allows to reduce to the situation where the energy 

r 2 



density in f)3. 15j) is bounded below. For the proposition below we define for all x' £ T? • 



u*(x') := £ \Af \Hf(x% (5.1) 

where 5f(x') is given by (|3.18p . We also recall that p £ defined in (|3.ip is the expected radius 
of droplets in a minimizing configuration in the blown up coordinates. 

We cover T| e by the balls of radius |ro whose centers are in ^Z 2 . We call this cover 
{U a }a and {x a } a the centers. We also introduce D a := B(x a , ^). 

Proposition 5.1. Let h' e satisfy (|3.16p . assume (|2.19p holds, and set 

f £ := \Vh'\ 2 + — \h'\ 2 -— \\np £ \u e . (5.2) 
js i e i 2 |lne|' el 2vr' P 1 v ; 

Then there exist £q > as in Proposition and constants c, C > depending only on 5 
and k such that for all e < e$, there exists a family of integers {n a } a and a density g £ on 
T 2 e with the following properties. 

- g £ is bounded below: 

g £ > -cln 2 (M e + 2) onT%. 

- For any a, 

n 2 a <C{g £ (D a ) + cln 2 (M £ + 2)). 

- For any Lipschitz function \ on ^% we have 

X(fe ~ 9e)dx' 



< CY,(v £ (U a ) + K + M £ ) In(n Q + M £ + 2)) ||VxlU<- (Do) . 

a 

(5.3) 
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Proof. The proof follows the method of [32] , involving a localization of the ball construction 
followed by energy displacement. Here we follow \35\ Proposition 4.9]. One key difference 
is the restriction to Ip which means we cover only those Vt^ satisfying Af > /3 as in 
Proposition (|4.2p . 

- Step 1: Localization of the ball construction. 

We use U a defined above as the cover on T? s . For each U a covering at least one droplet 
whose volume is greater or equal than /3 and for any r G (r(Bo), |ro) we construct disjoint 
balls £>" covering all e with i G 1/3 u a s using Proposition 14.21 Then choosing a small 
enough p G (r(^o) ; \fo) independent of e (to be specified below), we may extract from 
L) a Bp a disjoint family which covers Ui e / S f^ £ as follows: Denoting by C a connected 
component of U Q £>", we claim that there exists ao such that C C J7 ao . Indeed if x G C and 
letting A be a Lebesgue numbei0 of the covering of T| £ by {U a } a (it is easy to see that 
in our case jTq < A < ^f*o), there exists ao such that B(x, A) C ?7 Q0 . If C intersected the 
complement of U ao , there would exist a chain of balls connecting x to (U ao ) c , each of which 
would intersect U ao . Each of the balls in the chain would belong to some £>" with a' such 
that dist (U a i,U ao ) < 2p < ^rQ. Thus, calling k the universal maximum number of a"s 
such that dist (U a >, U ao ) < ^ro, the length of the chain is at most 2kp and thus A < 2kp. 
If we choose p < X/(2k), this is impossible and the claim is proved. Let us then choose 
p = A/(4/c). By the above, each C is included in some U a . 

We next obtain a disjoint cover of Ui e i„Q' iE from U a Bp. Let C be a connected component 
of U a Bp. By the discussion of the preceding paragraph, there exists an index ao such that 
C C U ao . We then remove from C all the balls which do not belong to Bp and still denote 
by Bp the obtained collection. We repeat this process for all the connected components 
and obtain a disjoint cover B p = U a Bp of Ui^i^' ie . Note that this procedure uniquely 
associates an a to a given B G B p , as well as to each for a given i G Ip by assigning 
to it the ball in S p that covers it, and then the a of this ball. We will use this repeatedly 
below. We also slightly abuse the notation by sometimes using Bp to denote the union of 
the balls in the family Bp. 

We now proceed to the energy displacement. 

- Step 2: Energy displacement in the balls. 

Note that by construction every ball in Bp is included in U a . From the last item of 
Proposition 14.21 applied to a ball B G Bp, if e is small enough then, for any Lipschitz 
non-negative x we have for some c > depending only on k and 5 and a universal C > 

> -Cv e (B)\\V X \\L~(B), 

1 A Lebesgue number of a covering of a compact set is a number A > such that every subset of diameter 
less than A is contained in some element of the covering. 
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where v e is defined by (|5.ip . Rewriting the above, recalling the definition (|3.ip and defining 
n a > 1 to be the number of droplets included in U a D B and satisfying A\ > f3, we have 

L * H 2 + iJ - i - c ) £ ^ + L ^ 

> -Cu e (B)\\Vx\\ L ~(B), 
where we set f a := ^° and define (recall that a implicitly depends on i E Ip) 

<*oo --=lE ( r ^f) ^ = I E i^i 2ln (t) ^ (5 - 4) 

Z7r i&Ip \ nape J ^ i£l VTW 

The quantity uj 6 in some sense measures the discrepancy between the droplets and 
balls of radius p e . We will thus naturally use M £ in (|3.13p to control it. Note also that it 
is only supported in the droplets, hence in the balls of B p . 
Applying Lemma 3.1 of |32] to 

we deduce the existence of a positive measure gs,e such that 

\\f B ,e-9B,ehip* <C^(B), (5.5) 

where Lip* denotes the dual norm to the space of Lipschitz functions and C > is universal. 
- Step 3: Energy displacement on annuli and definition of g £ . 

We define a set C a as follows: recall that p was assumed equal to A/(4fc), where A < \tq 
and k bounds the number of a"s such that dist (U a >, U a ) < |ro for any given a. Therefore 
the total radius of the balls in B p which are at distance less than r$ from U a is at most 
kp = jqTq. In particular, letting T a denote the set of t G ^p) such that the circle of 
center x a (where we recall x a is the center of U a ) and radius t does not intersect B p , we 
have \T a \ > -^r^. We let C a = {x \ \x — x a \ G T a } and recall that D a = B(x a , ^p). 
Let t G T a . Arguing exactly as in the proof of (|4.10p . we find that 

I \ Vh 'f d u 1 (x') + -^-[ \hXdx'>^(l--^-] 

JdB(x a ,t) M^^\JB( Xa ,t) - 2irt \ 4|lne|7 

with m £: t '■= f B ( x t ){^' £ {, x ')~ P e ) dx 1 . Arguing as in (|4.12p and using the fact that B(x a , ^ro) 
contains all the droplets with i G IpUai we nn d that we can take e sufficiently small 
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depending on k, and tq sufficiently small depending on k and 5 such that for all t £ T a 

[ \ Vh ' £ f d u 1 (x') + -!f- [ \h' e \ 2 dx'>-L( V Af). 

JdB(x„.t) 4|lne| J B(Xn , t] l 4vrtl J 



Integrating this over t G T a , using that \T a \ > t^o, we obtain that 



■ 



f |V^|W + -g- / \h>fd X >>c( Y, V) , (5-6) 

with c > depending only on rg, hence on k and 5. 

We now trivially extend the estimate in (|5.6p to all a's, including those U a that contain 
no droplets of size greater or equal than f3. The overlap number of the sets {C a } a , defined 
as the maximum number of sets to which a given x' 6 Tf e belongs is bounded above by 
the overlap number of the sets {D a } a , call it k' . Since the latter collection of balls covers 
the entire T^ e , we have k' > 1. Then, letting 



Be Bp 



and 

.2 



'-WK + S 1 ^^^ (in £-c)§-. (5.8) 



we have 



/;-E^s=(l^l 2 + ^KI 2 )i v ^ 

-W?(lv^l 2 + ^KI 2 )i C . + ^KlV 



and from (15.6 



^™ = w/ c H 2 + ^<K + i( ln £- c ) £ AT -MA.) 
>J Yl A i) -^ lnn ° E iifi 2 -w e p«)-c E i^i 2 ' ( 5 - 10 ) 
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for some C, c > depending only on k and 5. Now we combine the middle two terms, 
using the definition of oj £iCI in (|5.4|) . to obtain 

fa,e(Da) > J Y, ^) -^ lnf " E \ A I\ 2 ~ C E l^l 2 ' ( 5 - U ) 

The next step is to bound f a . We separate those with Af > 3 2 / 3 7T7 -1 and those 
with Af < 3 2 / 3 7T7 _1 . We denote (with s for "small" and 6 for "big") 

IKa = {i£ Ip,u a ■■ Af < S^ttt- 1 } , 
For the small droplets, we use the obvious bound 

Y i^i 1/2 ^ m ^ ( 5 - 12 ) 

with a universal c > 0, while for the large droplets we use that in view of the definition of 
M £ in (gJ5D we have 

Y |Af I 1 / 2 < C ^ Af <C'M £ , (5.13) 

P,a p,ck 

for some universal C,C > 0. We can now proceed to controlling f a . By (|3.ip and (|4.3p . 
for universally small e we have 

r a <C Y Pi, (5-14) 

for some universal C > 0. In view of (|3.13p . ()5. 13j) and (|5.12p . we deduce from Remark 
that for universally small e we have 

<C(M £ + cn as + C7'M e ) < C"{n as + M £ ) < C"(l + n as + M £ ), (5.15) 
where c,C,C ,C" > are universal. Therefore, (|5.1ip becomes 

2 / \ 2 



f a , e (D a )>c( Y A t] + c ( E A i) - Clnf - E \ A i\ 2 - c '" E i^i 2 ' 

>c0 2 < +t / aA -C'ln(C" / (l + n as +M £ ))[n as + £ AfY (5.16) 
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where C, C > are universal, c, C" , C" > depend only on k and 5, and C was chosen 
so that C\AI\ 2 < C\A\ + 1). 

We now claim that this implies that 




WO tt )>/< + d -C'"l^(M £ + 2), (5.17) 



where C" > depends only on k and 5. This is seen by minimization of the right-hand 
side, as we now detail. For the rest of the proof, all constants will depend only on k and 
8. For shortness, we will set X := ^2 iGl t Af. 

First assume n Qs = 0. Then (|5.16p can be rewritten 

fa,e(D«) > cX 2 - C"ln(C"(l + M e )))X, 

By minimization of the quadratic polynomial in the right-hand side, we easily see that an 
inequality of the form (|5.17p holds. Second, let us consider the case n Qs > 1. We may use 
the obvious inequality log(l + x + y) < log(l + x) + log(l + y) that holds for all x > and 
y > to bound from below 

C -f?nl s + °-X 2 - C'ln(C"(l + n as + M E ))(n as + X) > |/3 2 n^ + °-X 2 
- C(n as +X) - Cn as \n(n as + 1) - CX\n{n as + 1) - Cln(M £ + l)(n as + X). (5.18) 

It is clear that the first three negative terms on the right-hand side can be absorbed into 
the first two positive terms, at the expense of a possible additive constant, which yields 

C -fnl s + °-X 2 - C'HC"(n as + M £ ))(n Qs + X) 

> -^ 2 n 2 as + -X 2 - Cln(M £ + l)(n«. + X) - C. (5.19) 

Then by quadratic optimization the right hand side of (|5. 19[) is bounded below by — C In 2 (M £ + 
2) (after possibly changing the constant). Inserting this into fj5. 16f) . we obtain (15.17|) . 

We then apply [32l Lemma 3.2] over D a to f a>E + C"'\D a \- 1 ln 2 (M e + 2), where C" is 
the constant in the right-hand side of (|5.17j) . We then deduce the existence of a measure 
g a ^ e on T 2 e supported in D a such that g a ^ e > —C"'\D ol \~ 1 \n 2 (M £ + 2) and such that for 
every Lipschitz function \ 



X(fa,e ~ 9a,e) dx' 



Da 



< 2diam (D a )\\V X \\L^(D a )fa,e(Da) 

<Cln(n Qs +M £ + 2)||V X ||L-(D Q ) £ \M\ 2 , (5-20) 
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and we have used the observation that 



and (|5.15p to bound the negative part of / a£ . In particular, taking x = 1) we deduce, in 
view of (|5.17p . that 

2 

(5.22) 



9a,e(Da) = /«, e (A») > ^ 2 < s + \ ( £ Af\ - C7'"ln 2 (M £ + 2), 



from which it follows that 



g a ,e(Da) > c « + (#/ /3 , a ) 2 ) - C"'\n 2 {M £ + 2) > y n 2 a - C"'\n 2 (M £ + 2). (5.23) 



Recalling the positivity of gs e introduced in Step 2, we now let 



9e ■= ^2 9B,e + ^2 9a,e + [ f'e ~ ^ fa,e ) , 
BeBp a V a J 



(5.24) 



and observe that since f' e — ^2 a f a ,e is also non- negative by (|5.9p . and since ^2 a g a ,e is 
bounded below by — k'C"'\D a \~ l ln 2 (M e + 2), where, as before, k' is the overlap number 
of {D a } a , we have g £ > — cm (Afg + 2) for some c > depending only on k and 5, which 
proves the first item. The second item follows from (|5.23p . (|5.24p and the positiveness of 
g B , £ and (f' £ - ^ a / a , e ). 

- Step 4 : Proof of the last item. 

Using the definition of g £ in (|5,24p . for any Lipschitz x we have 

/ X9edx = V / X9B,edx' + V / x(ffa,£ ~ fa,e)dx' + / x/gtfcc'. 

Hence, in view of (|5.5p . (|5.7p and (|5.20p we obtain for some C > 



x(/e - ffeW 



B^B P 



/ X(ftB,e - fB,e)dx' 



T- 



X{9a,e - fa,e)dx' 



<C £ v £ (B)\\Vx\\l°°(b)+CJ2 Hn as +M £ + 2)\\V X \\L^ { D a ) £ l 1 ^- (5-25) 



Using that \A\ | 2 < C(t4| + 1) for a universal C > and (|5.13p . we have 

J2 \Af\ 2 <C(n as +M £ ). 



Since n Qs < n a , the third item follows from (|5.25p . 



□ 
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We now apply Proposition 15.11 to establish uniform bounds on M £ , which characterizes 
the deviation of the droplets from the optimal shape. 

Proposition 5.2. If (12.19|) holds, then M £ is bounded by a constant depending only on 
sup e>0 F £ [u £ ] , k, 5 and t. 

Proof. From the last item of Proposition 15.11 applied with x = \ together with the first 
item, we have 

/ f £ dx' = [ g £ dx' > -C|lne|ln 2 (M e + 2), 

with some C > depending only on k, 5 and £, while from (|2.19p . (|3.12|) and (|5.2p . we 
have 

C > £ 2 F £ [u £ ] > M £ + — ^— / f £ dx' + oil) > M £ - Cln 2 (M £ + 2) + oJl), 

l ln ^l Jt% 

for some C > depending only on sup e>0 F £ [u £ ], k, 5 and i. The claimed result easily 
follows. □ 

With the help of Proposition 15.21 an immediate consequence of Proposition 15.11 is the 
following conclusion. 

Corollary 5.3. There exists C > depending only on k, 5, £ and sup £>0 F £ [u £ ] such that 
tf 9e is as i n Proposition I5.il and (|2.19p holds, then g £ > —C. 

In the following, we also define the modified energy density g £ , in which we include 
back the positive terms of M £ and a half of p^l^el 2 that had been "kept aside" instead 
of being included in f £ : 

5 £ :=^ + ^^|/4| 2 + |lne|j J>? - y/^SJ+a £ A\h 

11 1 i A e >7r3 2 / 3 7- 1 



^ (A £ -7rf 2 ) 2 ~5 i + c 3 £ ^1 (5.26) 

<7r3 2/3^-l Af<8 J 



+ C 2 



where we recall r £ = ( )} ne \ ) and 5f is defined by (|3.18|) . These extra terms will be 



1/3 

used to control the shapes and sizes of the droplets as well as to control h' £ . We also point 
out that in view of (|5.2p . (|5.3|) and (|3.12p . we have 



" 2 F £ [u £ ]>-^— f g £ dx' + o £ (l). (5.27) 
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6 Convergence 

In this section we study the consequences of the hypothesis 

Vi? > 0, Cr := limsup / g £ (x + x £ )dx < +00, (6-1) 
e^O JK R 

where Kr = [— R, R] 2 and (x®) is such that x® + Kr C T| e . This corresponds to "good" 
blow up centers x £ , and will be satisfied for most of them. 

In order to obtain o e (l) estimates on the energetic cost of each droplet under this 
assumption, we need good quantitative estimates for the deviations of the shape of the 
droplets from balls of the same volume. A convenient quantity that can be used to char- 
acterize these deviations is the isoperimetric deficit, defined as (in two space dimensions) 

Wei 

D(Q' he ) := l *f_ - 1. (6.2) 

The isoperimetric deficit may be used to bound several types of geometric characteristics 
of that measure their deviations from balls. The quantitative isoperimetric inequality, 
which holds for any set of finite perimeter, may be used to estimate the measure of the 
symmetric difference between and a ball. More precisely, we have [T7] 

a(n' h£ ) < Cy^T), (6.3) 
where C > is a universal constant and aiVt^ ) is the Fraenkel asymmetry defined as 

a ( n i,e) : = mm |^/ 1 » ( 6 - 4 ) 

where A denotes the symmetric difference between the two sets, and the infimum is taken 
over balls B with \B\ = \£l' ie \. In the following, we will use the notation r\ and af for the 
radii and the centers of the balls that minimize a(£l' ie ), respectively. 

On the other hand, in two space dimensions the following inequality due originally to 
Bonnesen [6] (for a review, see [30]) is applicable to £l' i£ : 

Rt<ri(l + c^D(W l£ )). (6.5) 

Here Rf is the radius of the circumscribed circle of the measure theoretic interior of fi^ 
and c > is universal. Indeed, apply Bonnesen inequality to the saturation of Cl^ £ (i.e., the 
set with no holes) for each droplet. Then since the set £ is connected and, therefore, its 
saturation has, up to negligible sets, a Jordan boundary [4], Bonnesen inequality applies 
to it. 
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6.1 Main result 

We will obtain local lower bounds in terms of the renormalized energy for a finite number 
of Dirac masses in the manner of [5] : 

Definition 6.1. For any function x an d f> £ A m (cf. Definition \2.1\) , we denote 

W(<p,x) = lim [ \ [ x|Vv?| 2 dx + vrln?7^x(p) ) • (6-6) 



r?->0 I 2 7R2\ Up6A B(p,r7) 



We now state the main result of this section and postpone its proof to Section 
Throughout the section, we use the notation of Sec. To further simplify the notation, 
we periodically extend all the measures defined on T| e to the whole of R 2 , without relabeling 
them. We also periodically extend the ball constructions to the whole of R 2 . This allows 
us to set, without loss of generality, all x® = 0. 

Theorem 4. Under assumption (|2.19p . the following holds. 

1. Assume that for any R > we have 

\i-msxxp g £ {K R ) < +oo, (6.7) 

where Kr = [— R, R] 2 . Then, up to a subsequence, the measures fi' e , defined in 
(|3.17p . converge in (Co(R 2 ))* to a measure of the form v = 3 2 / 3 7r^ agA 5 a where A 
is a discrete subset o/R 2 , and {(p e } £ defined in (|2.17p converge weakly in W^(M. 2 ) 
for any p G (1,2) to (p which satisfies 



-A(^ = 2vr^5 a 



m in K , 

aGA 

in the distributional sense, with m = 3 _2 / 3 (5 — 5 C ). Moreover, for any sequence 
{f2j e>e } £ which remains in Kr, up to a subsequence, the following two alternatives 
hold: 

i. Either A3 < £^ and Pf < -SflL- as e -»• 0, 

ii. Or A\ s is bounded below by a positive constant as e — > 0, and 

A5 -> 3 2/3 vr and Pf -> 2 • 3 1/3 vr as e — > 0, 

with 

a(nl >£ )<— ^ ase^O, (6.8) 
me ' 
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for some Cr > independent of e. 
2. If we replace (|6.7|) by the stronger assumption 



lim sup g £ (K R ) <CR 2 , (6.9) 
where C > is independent of R, then we have for any p E (1,2), 

limsup (r— | / \Vip\ p dx J < +oo. (6.10) 



Moreover, for every family {xr}r>o defined in Definition \2.3\ we have 

a6A 



liminf / XRQedx > -—W(<p, X r) + V Xfl(a) + o(|^fl|). (6.11) 

s->o J R2 2 8 f-f 



Remark 6.2. W^e poini oui that it is included in Part 1 of Theorem^ that at most one 
droplet Q! ie s mft ^4j e . e bounded from below converges to a 6 A. Indeed otherwise in the first 
item we would have fj,' £ — > 3 2 / 3 irn a X^ a gA ^ a where n a > 1 is the number of non-vanishing 
droplets converging to the point a. 

Theorem 2] relies crucially on the following proposition which establishes bounds needed 
for compactness. Each of the bounds relies on (|6.7p . Throughout the rest of this section, 
all constants are assumed to implicitly depend on k, 5, I and sup e>0 -F £ [it e ]. 

Lemma 6.3. Let g e be as above, assume (|6.7p holds and denote Cr = limsup^g (^(-JCr) . 
Then for any R and e small enough depending on R we have 

]T n 2 a <C(C R+C + R 2 ), (6.12) 

a \u a CK R 

J2 A!<C(C R+C + R 2 ), (6.13) 

I XR(fe-9e)dx <C V K + l)ln(n a + 2) <C{C B+C + R 2 ), (6.14) 

Jk r i 

a \u a CK R+c \K R _ c 

where {xr} is as in Definition \2.S\ and n a = jflpjj a , with U a as in the proof of Proposition 
5.11 for some C > independent of e or R. Furthermore, for any p G (1, 2) there exists a 
C p > depending on p such that for any R > and e small enough 

[ \Vh' £ \Pdx < C P (C R+C + R 2 ). (6.15) 

JK R 
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Proof. First observe that the rescaled droplet volumes and perimeters A\ and Pf are 
bounded independently of e, as follows from Proposition 15.21 and the definition of M £ . 
Then, (|6.12p and ()6.13|) are a consequence of (|6.7p . the second item in Proposition 15.11 
together with the upper bound on M e . The first inequality appearing in (|6.14p follows 
from item 3 of Proposition 15.11 with the bound on M e , where we took into consideration 
that only those D a that are in the O(l) neighborhood of the support of | Vxr| contribute 
to the sum, along with the observation that the mass of v e (of (|5.ip ) is now controlled 
by n a (a consequence of the above fact that all droplet volumes are uniformly bounded). 
The second inequality in (|6.14p follows from (|6.12p . The bound (|6.15p is a consequence 
of Proposition 14.21 and follows as in [32] and [35]. We refer the reader to [32], Lemma 4.6 
or [35] Lemma 4.6 for the proof in a slightly simpler setting. □ 

6.2 Lower bound by the renormalized energy (Proof of Theorem 3]) 

We start by proving the first assertions of the theorem. 



- Step 1: All limit droplets have optimal sizes. ^From (|5.26p . (|6.7p and Corollary I5.3[ for 
all e sufficiently small depending on R we have 



JKr V i V 7 A|>3 2 /37r 7 -i 

£ (Af - vflfll + c 3 J2 m) dx < ^ , (6.16) 

<7r3 2/3 7 -i A!<0 J 



/3<Af<T. 

where we recall that all the terms in the sums are nonnegative. It then easily follows that 
for all i £ Ik r the droplets with Af > 3 2 / 3 7T7 _1 do not exist when e is small enough 
depending on R, and those with A\ < (3 satisfy A\ = Cr\ lne| _1 and Pf < Cr\ lnep 1 ^ 2 , 
for some Cr > independent of e. This establishes item (i) of Part 1 of the theorem. 

It remains to treat the case of Af £ [j3, 3 2 / 3 7T7 _1 ] when e is small enough. It follows 
from (flTToT) that 

D(ftJ < , (6.17) 

for some Cr > independent of e, and since f e = 3 1 / 3 + o e (l), for all these droplets (or 
equivalently for all droplets with Af > (3) we must have 

^^3 2/3 vr and Pf -»• 2 • 3 1/3 vr as e -> 0. (6.18) 

Using dS3J, §61$ easily follows from (1638]) and (|B1B|> . 
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- Step 2: Convergence results. 

^,From boundedness of Af, (|6.13p and (|6.16|) we know that #Ip,K R an d I^' £ (Kr) are both 
bounded independently of e as e — > 0. We easily deduce from this, the previous step and 
the definition of fi' £ that up to extraction, fj,' £ converges in each Kr to at most finitely 
many point masses which are integer multiples of 3 2//3 7r and, hence, to a measure of the 
form v = 3 2 / 3 7r XlaeA da$a, where d a G N and A is a discrete set in the whole of M 2 . In 
view of (|6.15p . we also have h £ h G W^(M. 2 ) as e — > 0, up to extraction (recall that we 
work with equivalence classes from (|2.10j) ) . Finally, from the definition of g e in (|5.26j) and 
the bound (16.71) we deduce that 



K 2 

In el 



/ \K\ 2 < c R 

JKr 



from which it follows that | lne| _1 /i^. tends to in L 2 oc (IR 2 ) as e — > 0. Passing to the limit 
in the sense of distributions in (|3.16p . we then deduce from the above convergences that 
we must have 

- Ah = 3 2/3 vr d a$a - p on M 2 . (6.19) 

aeA 

We will show below that d a = 1 for every a G A, and when this is done, this will complete 
the proof of the first item after recalling (p £ = 2 • 3~ 2 ^ 3 h' £ and m = 2 • 3~ 2//3 /L 

- Step 3: There is only one droplet converging to any limit point a. 

In order to prove this statement, we examine lower bounds for the energy. Fix R > 1 such 
that 8Kr n A = and consider a G A n Kr. From Step 1, (|2.2p and Lemma 14.11 for 
any n G (0, ^) such that r\ < \ miiibeAriK R \{a} \ a ~ b\ and for all r < n, all the droplets 
converging to a are covered by B(a, r), and B(a, rj) contains no other droplets with Af > f3, 
for e small enough. There are d a > 1 droplets in B(a,r) such that A\ — > 3 2//3 7r as e — > 0, 
let us relabel them as fi^ e , . . . , Q' d £ . 

Let U = B(a, rj). Arguing as in the proof of the first item of Proposition 14.21 by (16.181) . 
we may construct a collection Bq of disjoint closed balls covering Uie/^ v e an< ^ satisfying 

r(Bo) < Cd aPe < n, (6.20) 

for some universal C > 0, provided e is small enough, and a collection of disjoint balls B r 
covering Bq of total radius r G [r(Bo),rj\. Choosing r = n 3 , which is always possible for 
small enough e, it is clear that B^s consists of only a single ball contained in B(a, |?? 3 ) for 
e small enough. Applying the second item of Proposition 14.21 to that ball, we then obtain 
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Therefore, we have 



4 » (l v '<l 2 + * s i ('" 1W) - ^) fe> ") t |if " <6 ' 22) 

On the other hand, we can estimate the contribution of the remaining part of B(a,rj) as 
/ X R\Vti e \ 2 dx' + -^- I XR\K\ 2 dx' 

JB{a,r))\B v 3 4 I ln£ \ JB{a,ri) 

> ( min xr] ( [ \VK\ 2 dx' + / \h'f dx') 

\B(a,r)) J yj B(a, V )\B(a,2r)3) 4 I ln e \ J B(a, V ) J 



> I min xii I / 1/ |V^| 2 d«V) + 7ri— T / |^| 2 ds' dr B . (6.23) 

\B(a,r)) J J2ri 3 \JdB{a,r B ) 4 I m £ l JB{a,r B ) J 

Arguing as in (|4.12p and using the fact that rj < ^, we obtain 

X(a,7,)\B n3 e ' ' 4 l ln£ l 



K 2 

XR \Vti £ \ 2 dx> + — r / XR\K\ 2 dx' 

B(a, V ) 

1 

> — I min yr 

~ 27T \B(a,r,) A 



) ln ^(f>*) C 1 -^ ( 6 - 24 ) 



where C > is independent of ry and e, for small enough e. 

We will now use crucially the fact shown in Step 1 that all Af > (3 approach the same 
limit as e — > 0. We begin by adding (|6.2ip and (|6.24p and subtracting lnp e | X)j=i 1^1 | 2 Xr 
from both sides. With the help of (|6.20p we can cancel out the leading order 0{\ ln p e \) term 
in the right-hand side of the obtained inequality. Replacing Af and Af with 3 2 / 3 7r + o e (l) 
in the remaining terms and using the fact that ming( a ^ XR > XR( a ) ~ ^v\\^Xr\\oo on 
B(a,r]), we then find 



f ( i. — , . / iO ft i , / iO 

J B{a,r]) 



XR ( IV/l'l 2 + — -\h'\ 2 - —I lnpJi^ ) dx' 

el 2 hie 1 £l 2vr' 



>— X*(a)(^ln^ + d a ln-LJ-C', (6.25) 

where C > is independent of e or r\. 

Now, adding up the contributions of all a 6 A n and recalling the definition of f £ 
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in (|5.2p . we conclude that on the considered sequence 
limsup/ XRfedx' > limsup V] / XFtfedx' 

e^O J Kr e^O aeAnKR JB(a, V ) 

>^|l nr? | £ ( 2c £ - 3d ) X ii(a) - C, (6.26) 

for some C > is independent of e or 77. In particular, since xn( a ) > for all a G A n Kr, 
the right-hand side of (|6.26p goes to plus infinity as rj — > 0, unless all d a = 1. But by 
the estimate (|6.14p of Proposition 16.31 Corollary 15.31 and our assumption in (|6.7p together 
with (|5.26p . the left-hand side of (|6.26p is bounded independently of rj, which yields the 
conclusion. 

- Step 4 : Energy of each droplet. Now that we know that for each aj 6 An Kr there exists 
exactly one droplet e such that a\ — > ai and A\ — > 3 2//3 7r, we can extract more precisely 
the part of energy that concentrates in a small ball around each such droplet. Let B^ be a 
ball that minimizes Fraenkel asymmetry defined in (|6.4p . i.e., let Bi = B(af,rf), and let B 
be a ball of radius re centered at a\ . Arguing as in (|4.12p in the proof of the second item 
of Proposition 14.21 we can write 

\Vh | W(x') + — - / h' £ 2 dx' > '^—^ -(I -erf). (6.27) 

dB 4|lne|y B 2nr B 

Observe that by the definition of Fraenkel asymmetry we have |^ e ni?| > \B\ — ia(^ £ )|-Bj| 
for all re < rf. Hence, denoting by rf the smallest value of rs for which the right-hand 
side of this inequality is non-negative and integrating from ff to rf , we find 

(\vK\ 2 + 7] ^-\K\ 2 )dx' 

Bi V 4|lne|' J 

> ^(1 + 0e (l))e-V3| lne |-2/3 r i r -l (r |_| ff |2 )2drB _ (62g) 
2 Jr e 

Since by (|6.8p and f|6. 18j) we have rf/rf — > and e~ 1//3 | lne|~ 1//6 r[ — > 3 1//3 as e — > 0, after 
an elementary computation we find 



X,,(i v ^' 2+ ^"<' 2 ) di ' a ^ +oc(l) - (6 - 29) 

On the other hand, by (|6.5p and (|6.17p it is possible to choose a collection £> C B(ai,rj), 
actually consisting of only a single ball B(df,Rf) circumscribing £l' ie , so that 

r(Bb) = #f < ^ (l + Cfl| lnel" 1 / 2 ) = Pe + o e (p e ). (6.30) 
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The corresponding ball construction B r of the first item of Proposition 14.21 with U = 
B(ai,rj) and r\ as in Step 3 of the proof (again, just a single ball B(af,r)), exists and is 
contained in U for all r € [r(Bo),r)'], for any rf € (r(Bo),rj), provided e is sufficiently small 
depending on rf . In view of the fact that for small enough rf and small enough e depending 
on rf we have xr( x ) > Xii(^f) ~~ c \ x ~ ®i I > 0, with c > independent of e, 77' or R, we 
obtain that 

f XRNK\ 2dx ' + JTT 1/ Xfil^el 2 ^' 

J B(al,r,')\B 4 I™ £ I ^(af^) 

> f (Xfl(5f ) - cr) f / |V/4| WO*)") * + / l^l 2 ^' 

Jr(Bo) \J,9B r / 8|ine| JB(af,r]') 

> f V (Xdaf) - cr) ( I \Vti £ \ 2 dU l (x) + [ \h' £ \ 2 dx') dr 

Jr(Bo) \JdBr °V I [n£ \ JB(af, v ') J 

>T IXiffi) -<*)(( \Vh' e \ W(x) + [ \h' £ \ 2 dx')dr 

Jr(Bo) \JdBr 4 I me l JBr J 

1 n' dr 

> ^\Af\ 2 / (XR(a!) - cr)(l - Cr)-, (6.31) 

for rf and e sufficiently small, arguing as in (|4.12p in the proof of Proposition 14.21 and 
taking into account Remark 14. 41 in deducing the last line. Performing integration in (I6.31|) 
and using (|6.30p . we then conclude 

/ XR \ Vh i £ f dx > + ^L— ( X R\K\ 2 dx' 

J B(aS,rf)\Bo 4 I m e UB(af ,77') 

>i-|if| 2 Xfl (af)ln(^)-C ? /, (6.32) 

for e sufficiently small. 

- Step 5: Convergence. Using the fact, seen in Step 2, that h' e h in W^iM 2 ), we have, 
by lower semi-continuity, 

liminf / XR\^h' £ \ 2 dx > / XR|V/i| 2 £te'. (6.33) 
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On the other hand, in view of xr(5|) = Xr + 0(p e ) by (|6.30p . from (|6.32p we obtain 

liminf f XRNKfdx' + [ XR ( J]^\K\ 2 ~ ^p £ \v £ ) dx' 

JB(af,ri)\B J B(af,r,) \ 4 l me l 27r J 

> liminf / X R\Vti £ \ 2 dx' + / XR ( Tn~^\K\ 2 ~ ^p £ \v £ ) dx' 

JB(af,ri')\Bo JB(a\,rf) V 4 l lne l 27T J 

3 4 / 3 vr 

> —XRia^'-Cr,', (6.34) 

where we also used that Xr ~^ XR{ a ) as e — > 0. 

We now convert the estimate in (|6.29|) to one over Bq and involving xr as well. Ob- 
serving that Q' i£ C Bq and that xr( x> ) > xf ~ 4/0 e ||Vxi?||oo for all x' S £l' ie and e small 
enough by (|6.30p . from f|6.29[) and (|3.ip we obtain 

liminf / XR ( iV/i'J 2 + j^IKA dx' > ^XrK), (6.35) 
J Bo V 4|lne| J 8 

where we used the fact that by (|3.2|) . (|3.12p and (|4.5p the integral in the left-hand side 
of (|6.29p may be bounded by C|lne|, for some C > independent of e and R. Adding 
up (|6.33p with (|6.34p and (|6.35p summed over all a, £ Kr, in view of the arbitrariness of 
rj' < rj we then obtain 

liminf / X R (\Vh' £ \ 2 + ^^Kl 2 - i lnp £ \A dx' 

s^O J R 2 \ 2|hl£| 27T / 

[ XR \ Vh f dx ' + tllJ2xR(a)(]n V +])-C V . (6.36) 

JB*\U aeA B(a,Ti) 2 t^A ^ A ' 



> 



Letting now r\ — > in (|6.36p . and recalling that ip = 2 • 3 2 / 3 h and that the definition of 
W(ip, x) is given by Definition 16 .14 we obtain 



liminf f XR (\Vh' £ \ 2 + -^-\h' £ \ 2 -^\lnp £ \Adx' 
e->o J R 2 \ 2|me| 2tt J 



q4/3 

> —W(^xr) + —-Y.*^- ( 6 - 3? ) 



2 

agA 

iProm (I6TT41) we may replace f e = \S/h' £ \ 2 + 2TEJ\\K\ 2 ~ h\ ^p e \v e by g e in (15371) with 



an additional error term: 



r 3 4 / s 7r 
liminf / x^dx' > — W(<p, X r) + Y, Xr(*) ~ cA(R), (6.38) 

e^O J R2 2 8 ^ 
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where 

A(R) = limsup (n Q + l)ln(n Q + 2), 

£->0 

for some c, C > independent of i2. Under hypothesis (|6.9p , from (|6.12p we have 

limsup V n 2 a <CR 2 , 

a \U a CK R 

and thus, using Holder inequality and bounding the number of a's involved in the sum by 
CR we find 



A(R) < C limsup (n^ 2 + 1) 



< C'i? 1 / 4 limsup 



\ a \u a 



3/4 

n 2 I + CR < C"R 7 /\ 



^R+C 



for some C, C , C" > independent of R. Hence 



limsup limsup — = 0, 
which together with ()6.38p and the fact that g £ > g £ establishes (|6.1ip . □ 

6.3 Local to Global bounds via the Ergodic Theorem: proof of Theorem 
[H item i. 

The proof follows the procedure outlined in [35] . We refer the reader to Sections 4 and 6 
of [35] for the proof adapted to the case of the magnetic Ginzburg-Landau energy, which 
is essentially identical to the present one, with some simplifications due to the fact that 
we work on the torus. As in [35j . we say that n £ .A/fo(]R 2 ), if the measure d^i + Cdx is a 
positive locally bounded measure on K 2 , where C is the constant appearing in Corollary 
15.31 The measures dg £ and the functions tp £ will be alternatively seen as functions on T| £ or 
as periodically extended to the whole of M 2 , which will be clear from the context. We let x 
be a smooth non-negative function on M 2 with support in B(0, 1) and with J R2 x( x )dx = 1. 
We set X = W^(R 2 ) x M (R 2 ), and define for every x = (ip,g) E X the following 
functional 

f(x) := 2 / x (y)dg(y). (6.39) 
We note that from (|5.27j) we have for e sufficiently small 

F e [u £ ] + o £ (l) > [ dg £ = l f(9 x x £ )d\, (6.40) 
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where x £ := (if 6 , g e ), 9\ denotes the translation operator by A G ffi 2 , i.e., 9\f(x) := f(x+X), 
and -f stands for the average. Here the last equality follows by an application of Fubini's 
theorem and the fact that L 2 x( x )dx = 1. 

It can be easily shown as in |35] that f £ = f satisfies the coercivity and r-liminf proper- 
ties required for the application of Theorem 3 in [35] on sequences consisting of x £ = (tp e , g £ ) 
obtained from (u £ ) obeying (12.19|) . This is done by starting with a sequence {x £ } £ in X 
such that 

limsup / f(6» A x £ )dA < +oo, (6.41) 

e^O JK R 

for every R > 0, which implies that the integral is finite whenever e is small enough. 
Consequently f £ (# A x £ ) < +oo for almost every A € K^. Applying Fubini's theorem again, 
(I6.4ip becomes 

limsup / XR(y)dg~e(y) < +oo, 

where \R = X * 1k r , and "*" denotes convolution. Then since xr = 1 in Kr-\ and g £ is 
bounded below by a constant, the assumption (|6.7p in Part 1 of Theorem[5]is satisfied, and 
we deduce from that theorem that ip £ and g e converge, upon extraction of a subsequence, 
weakly in W^(M?) and weakly in the sense of measures, respectively. Furthermore, if x £ — > 
x = (</?, g) on this subsequence, we have 2 f R2 x(y)dg £ (y) = f (x £ ) ->■ f (x) = 2 J R2 x{y)dg{y). 

We may then apply Theorem 3 of [35] to f on T| e and conclude that the measure {P e } e 
defined as the push-forward of the normalized uniform measure on T| e by 

A ^ (6 X (p £ ,6 x g £ ), 

converges to a translation-invariant probability measure P on X with 

liminf F e [u £ ] > J f(x)dP(x) = J f*(x)dP(x), (6.42) 

where 

f*(<p,g)= lim / f(0 A x)dA= lim ( ^- [ X R(v)dg(y)) , (6.43) 

R^c°Jk r R^+co \\K R \ J R 2 J 

provided that x is in the support of P. 

The next step is to show that for P-a.e. x we have ip € A m with m = 3 -2//3 ((5 — 5 C ), and 
f* can be computed. By [351 Remark 1.6], we have that for P-a.e x, there exists a sequence 
{A e } £ such that x £ = (0\ s (p £ , 9\ e g £ ) converges to x in X. In addition, from (|6.42p ~ (|6.43p . 
for P-a.e. x, we have 

lim f {(0\x)d\ < +oo, 

R-,+ooJ Kr 
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for P-almost every x. Using Fubini's theorem again, together with the definition of f , we 
then find 



Therefore, since 



lim ( hfi I xn{y)dg(y) ) < +00. 

^+00 \ \K R \ J R 2 J 



xn{y)dg £ {y) ->■ / xn(y)dg(y) as e ->■ 0, 

JR2 



(6.44) 



a bound of the type (|6.9[) holds, and the results of Part 2 of Theorem U] hold for x E . In 
particular, we find that 

-Av9 = 2vr^5 a -m, (6.45) 

aeA 

with m = 3~ 2 / 3 ((5 — 5 C ), and that 

{*(<p,g) = lim lim / xji&tfc) > 3*/ 3 W(<p) + ^m. (6.46) 

fl^oo \\K R \ e->0j R2 ) 8 

The result in ()6.46p follows from the definition of f*, (|6.44p . (|6.1ip . the definition of W, 
provided we can show that 

Hm T±-Yx R (a)= lim ^4 = ™. (6.47) 

The latter can be obtained from (I6.15p . exactly as in Lemma 4.11 of [35], so we omit the 
proof. Note that with (|6.45p . it proves that <p £ A m , and we thus have the claimed result. 
Combining (j6.42j) and ()6.46p . we obtain 



£->0 



liminfF e [u e ] > I ^>' 6 W(ip) + ?^-(S-8 c ) \ dP(tp,g). 



Letting now P e and P be the first marginals of P £ and P respectively, this proves (12.201) 
and the fact that P-almost every ip is in A m with m = 3~ 2 / 3 (5 — 5 C ). □ 



7 Upper bound construction: proof of Part ii) of Theorem [T] 

We follow closely the construction performed for the magnetic Ginzburg-Landau energy 
in [35], but our situation is somewhat simpler, since we work on a torus (instead of a domain 
bounded by a free boundary). The construction given in [35] relies on a result stated as 
Corollary 4.5 in [35], which we repeat below with slight modifications to adapt it to our 
setting. These results imply, in particular, that the minimum of W may be approximated 
by sequences of periodic configurations of larger and larger period. Below for any discrete 
set of points A, |A| will denote its cardinal. 
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Proposition 7.1 (Corollary 4.5 in |35j). Let p G (1,2) and let P be a probability measure 
on W,'P (M?) which is invariant under the action of translations and concentrated on Ai- 
Let Q be the push- forward of P under —A. Then there exists a sequence R — >■ oo with 
R 2 G 2-7rN and a sequence {&r}r of 2R-periodic vector fields such that: 

- There exists a finite subset Ar of the interior of Kr such that 

( -div b R = 2vr Y. a e\ R S <* ~ 1 in K R 
\ bR ■ v = on OKr. 

- Letting Qr be the probability measure on Wj~ ' P (R 2 ), which is defined as the image 
of the normalized Lebesgue measure on Kr by x i-> —div 6_r(x + •), we have Qr — > Q 
weakly as R — > oo. 

- lim sup — \— lim \ \ [ \b R \ 2 dx + 7r|A R | hit] ) < / W{tp)dP{tp). 

R^oo \-Kr\ \2 JK R \U a eA R B(a,ri) J J 

Remark 7.2. We would like to make the following observations concerning the vector field 
bR constructed in Proposition \7.1\ 

1. By construction, the vector fields bR has no distributional divergence concentrating on 
8Kr and its translated copied since bR - v is continuous across 8Kr . However, bR ■ r 
may not be, and this may create a singular part of the distributional curl&#. This 
is the difficulty that prevents us from stating the convergence result for P directly in 
Theorem^ Part ii). 

2. We also note that an inspection of the construction in f35^ shows that bR is curl-free 
in a neighborhood of each point a G Ar and that curl&R belongs to W^' P (M?) for 
p < oo. 



7.1 Definition of the test configuration 

We take R the sequence given by Proposition IT. 11 The first thing to do is to change the 
density 1 into a suitably chosen density m Ej R, in order to ensure the compatibility of the 
functions with the torus volume. Recalling that \x £ > for 5 > S c and e small enough, we 
set 



m e ,R 



AR 2 



e|2 



2Rf e 



(7.1) 



where, as usual, [^J denotes the integer part of a x. We note for later that 



2fL £ 



CR , . 
<-=- = *(!)■ 



(7.2) 
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Recalling also that r £ = 3 1 / 3 + O ( ^rj^r ) and fjf — \x = O ( *" \^.f ) , we deduce that 



m £t R — > m, where m := 2 • 3 _2 / 3 /Z, as e — )• 0, for each i?. In particular, m £t R is bounded 
above and below by constants independent of e and R. The choice of m £) R ensures that we 
can split the torus into an integer number of translates of the square Kri with R' := -J= = , 

each of which containing an identical configuration of ^p- points. 

Let P € V be given as in the assumption of Part 2 of Theorem [TJ i.e., let P be a 
probability measure concentrated on A m . Letting P be the push-forward of P by ip \- > 
<-p{-j=)-i it is clear that P is concentrated on Ai, and by the change of scales formula (I2.16|) 
we have 



/ 



W((p) dP{ip) 



1 



m 



W(ip) dP(ip) + - lam. 



(7.3) 



We may then apply Proposition 17. II to P. It yields a vector field bR. We may then rescale 
it by setting 

b e ,R( x ) = V m £,RbR(.V m £,R x )- 



We note that b £t R is a well-defined periodic vector-field on Tj s because - ^^' R is an integer. 
This new vector field satisfies 



div b £ , R = 2ir ^ 5 a - m £)R in T 



(7.4) 



for some set of points that we denote A e> R, and 



1 



lim 

\K R \ v ^o 



V 



K R \U aFA - R B(a,r)) 



|6 £ijR J 2 (ix + vr|A e>jR n i^/^j-^l ln(r7 x /m £j ii) 



< 



J W(<p) dP(ip) + o R {\) as i? -> oo. 
Using (|7.3p and |A £i _r n K R j — | = this can be rewritten as 



^4 lim 
\Kr\ v^o 



1 

i * ./if g \Ua6A. R B{a,rfi 



Ifc^Pdx + TrlA^n^/^-^llnT? I -f— i^ln— ^ 



m 



< 



m 



W(<p) dP(cp)+o R (l) 



But we saw that m £: R — > m as e — > hence In (— £p) — > 0. Therefore, recalling the 
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definition of R' we have 



\^l l[m n\l I \b £ ,R\ 2 dx + 7r\A £iR nK R/vmr ^\\nri) 

\K R ,\ V ^0 ^2 JK R ,\U aeAeR B(a, V ) V ■ J 

< J W(<p)dP(<p)+o R (l) + o E (l). (7.5) 



It thus follows that 
1 



lim 



■ \ I Ift^lW + vrlA^llnT?) < f W(<p)dP(<p)+o R (l) + o e (l). 

{gey r^ y 2 hj e \U aeA£ R B(a, v ) J J 

(7.6) 

Note that A £jR is a dilation by the factor l/^/m e R , uniformly bounded above and 
below, of the set of points A R , hence the minimal distance between the points in A £ ^ R is 
bounded below by a constant which may depend on R but does not depend on e. For the 
same reason, estimates on b Rt£ are uniform with respect to e. 

In addition, we have that Q £)R , the push-forward of the normalized Lebesgue measure 
on T^ £ by x i y — div b £jR (x + -) converges to Q, the push-forward of P by —A, as e — > and 
R co. The final step is to replace the Dirac masses appearing above by their non-singular 
approximations: 

~5 a :=^^> r' £ :=e^\\ne\^f £) (7.7) 

where f £ was defined in (|3.ip . Note also that in view of the discussion of Section [3] it is 
crucial to use droplets with the corrected radius e l ^ \ In instead of its leading order 

value p £ = S^V^Ilnel 176 . 

Once the set A e>R has been defined, the definition of the test function u £ G A follows: 
it suffices to take 

u e (x) = -1 + 2 XB(a,r' e ) (x\ In e| 1/2 ) , 

which means (after blow up) that all droplets are round of identical radii r' £ and centered 
at the points of A £)R . We now need to compute -F £ [it £ ] and check that all the desired 
properties are satisfied. This is done by working with the associated function h' £ defined 
in (|3.16p . i.e. the solution in Tf £ to 

-Ati £ + -?—ti £ = Trf 2 £ (7.8) 

aeA e n 



obtained from (|3.16p by explicitly setting all A. 



- 7rr? 
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7.2 Reduction to auxiliary functions 

Let us introduce cp e , which is the solution with mean zero of 

- A(j) £ = 2ir ^ 5 a - m e! R in T% , 

where m £) R is as in (|7.ip . and f £ the solution with mean zero of 

- A/ e = 2tt 5 a - m £>R in T% . 



(7.9) 



(7.10) 



We note that f e is a rescaling by the factor m £) R — > m of a function independent of e, so 
all estimates on f £ can be made uniform with respect to e. 



Lemma 7.3. Let h' £ and <f> £ be as above. We have as e — > 



I \h' e \ 2 dx' <C R \\ne\ 



and for any 1 < q < oo 



v ( K - ^ £ 



< c R , 



Li(T%) 



(7.11) 



(7.12) 



for some constant Cr a > independent of e. 

Proof. Since A £ r is 2i?'-periodic, h' £ is too, and thus 



\h'\ 2 dx' = e\ 



\ne\f \h' e \ 2 dx' <C R \\ne\. 
Jk rI 



For the second assertion, let 

h £ (x) = h! e (x\f\ lne|) <j> e (x) = (/> e (xy/\ lne|) 
be the rescalings of h' £ and 4> £ onto the torus T|. Rescaling (|7.1ip gives 

IIMl 2 (t|) < Cr- 

Furthermore, the function w £ := h £ — \f 2 <\> £ is easily seen to solve 



(7.13) 



-Aw £ 



-K 2 ( h £ — 4- h F dx I 

V M J 



in T?. 



But from elliptic regularity, Cauchy-Schwarz inequality and (|7.13|) . we must have 



Ww £ 



< c 



f 



<c R . 



L2(T|) 



which yields (|7TT2j) . 



□ 
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The next lemma consists in comparing (j) £ and f £ . 
Lemma 7.4. We have 

Proof. We observe that f £ and <j) £ are both 2i?'-periodic. We may thus write 
<j> e (x)-f(x) = 2vr / G 2RI {x-y) V d(«5 a - 5 a )(y), 

where G*2_r' is the zero mean Green's function for the Laplace's operator on the square 
torus of size 2R' with periodic boundary conditions, i.e. the solution to 

-AG 2R/ =5 --^- 1 mT 2 2Rf (7.14) 

I 1 2R' I 

which we may be split as G2R'{x) = — log |x| + SW^r) with S2R' a smooth function. 
By Newton's theorem (or equivalently by the mean value theorem for harmonic functions 
applied to the function log | • | away from the origin) , the contribution due to the logarithmic 
part is zero outside of U ag A e R B(a,r' e ). Differentiating the above we may thus write that 
for all x g U a eA e R B(a,r' s ), 

V(0 e -/)(*) = 2tt / VSW(z-y) Y, d(6 a -5 a )(y). (7.15) 

/t2 — 

2R' a6A £ii? 

Using the C 2 character of S 2R ' we deduce that 

l|V(/ £ - ^)ll LTC (T| e \u 8 B(a,ri)) ^ c R'\Kr n K R ,\r' E 

and the result follows in view of (|7.7p . □ 

The next step involves a comparison of the energy of <f) e and that of b £t R and leads to 
the following conclusion. 

Lemma 7.5. Given A Et R as constructed above, and h' £ the solution to (17. 8ft . we have 



[ ^|V/4|W + 7r|A £iiJ |lnJ < f W(tp)dP(<p)+o e (l)+o R (l). 

Proof. In view of Lemmas 17.31 and 17. 4\ it suffices to show the corresponding result for 
f T 2 \i I B(an) \ fe\ 2 dx' instead of the one for h! e . From ()7.10|) and (|7.4p . we have 

div (6 £i _r — V/ e ) = hence by Poincare's lemma we may write V/ s = b eR + V -1 ^. We 
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note that — A£ e = curlb £j R, which is in W^' p for any p < +00 as mentioned in Remark 
17.21 By elliptic regularity we find that V£ e S L^ oc (IR 2 ) for all 1 < p < +00, uniformly with 
respect to e. We may thus write 

1 



/ T.\b e ,R\ 2dx ' 



= [ f^|V/ £ | 2 + l\V£ E \ 2 -Vf E -V L t e )dx', (7.16) 

^\UeA s ,/M) V 2 2 J 

where V/ e - V -1 ^ makes sense in the duality V£ e € L p , p > 2, V/ e G L q , q < 2. In addition, 
by the same duality, we have for any a £ A. e ,Ri 



lim / V/ £ • V^£ e = 

uniformly with respect to e. Therefore, we may extend the domain of integration in the 
last integral in (|7.16p to the whole of TT| e at the expense of an error 0^(1) multiplied by the 
number of points, and obtain 



i £ ,R\ 2 dx' 



[ \\Vfe?dx> < I \\b £ 

+ / V/ e -V- L £dx' + a fJ (|lne|). (7.17) 

Noting that the last integral on the right-hand side vanishes by Stokes' theorem (and by 
approximating V/ e and V -1 ^ by smooth functions), adding 7r|A £ ^| In?] to both sides, and 
combining with (|7,6p we obtain the result. □ 

In view of (I7,4j) and (I7.9P we have that — div b £ .R + Acp £ = 2ir X^aeA. R (^o, — $a) — > in 

W lo ^ ,p (M. 2 ), so we deduce, since the push-forward of the normalized Lebesgue measure on 
Tf s by x 1 ^ —div b £j ji(x + -) converges to Q, that the push-forward of it by x 1— > — Aip 6 (x + -) 
also converges to Q. Thus, part ii) of Theorem[T]is established modulo (|2.2ip . which remains 
to be proved. 

7.3 Calculating the energy 

We begin by calculating the exact amount of energy contained in a ball of radius 77. 
Lemma 7.6. Let h' e be as above. Then we have for any a £ A. £i r, 

/ \Vh/fdx' = — - + o e (l) (7.18) 

JB{a,r' £ ) ° 
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and 



L 



\Vh' e \ z dx' < -f e 4 ln^+o £ (l) + 0?? (l). (7.19) 

B(a,rj)\B(a,r' e ) * Pe 



Proof. In view of (|7.12p applied with q > 2 and using Holder's inequality, we have that for 
all a £ A EtR , 

r -2 2 

/ V(h' £ -^ £ ) dx'<o v {\). (7.20) 

JB(a,rj) z 

Thus it suffices to compute the corresponding integrals for <j> £ . Using again the 2R'- 
periodicity of (fi £ , we may write, with the same notation as in the proof of Lemma 17.41 



4>e{x) 



I G 2 r> (x - y) I 27r ^2 5 a (y) - m £ ,R.\ dy. 

\ aeA e , R J 



Since the distances between the points in A £< r are bounded below independently of e, and 
the number of points is bounded as well, we may write 4> e in B(a, rf) as 



cj) £ (x) = i/j £ (x) - / ln\x - y\ S a (y) dy (7.21) 



T 2 



where tp e ( x ) is smooth and its derivative is bounded independently of e (but depending on 
R). 

Thus the contribution of t/) £ to the integrals J^r arj ^ |V0 e | 2 is o v {\), and its contribution 
to /g( ar /) |Vi^> e | 2 is o £ (l). There remains to compute the contribution of the logarithmic 
term in (|7.2ip . But this is almost exactly the same computation as in (|6.27p - (|6.29p . and 
with (|7.20p it yields (17. 18f) . while it yields as well that 

/ \V(f) £ \ 2 dx' < 2vrln4 + o v {l). (7.22) 

J B(a, V )\B(a,r' e ) r s 

Now 

rl 1 f\lne\\ 1 ^ 3 ( ^/ln|lne|^ 

1 + 



p £ 3V3 VI In p e | 7 V V line] 

Consequently In — = o e (l), and so we may replace r' £ with p £ at an extra cost of o £ (l) in 
(17321) . and the result follows with (f7T20l) . □ 

We can now combine all the previous results to compute the energy of the test-function 
u £ . By following the lower bounds of Proposition 13. 1\ it is easy to see that in our case (all 
the droplets being balls of radius r' £ ) all the inequalities in that proof become equalities, 
and thus recalling (13. If) : 



K 



2 



F E [u £ ] = T7^ 2 / [\Vh> £ \ 2 + -—\ti £ \ 2 \dx' + Kf A £ \k £ , R \\np £ \+o £ {l 



I Jt* \ £l |ln e l' £l I 
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with the help of Lemma 17.61 we have for every R 

i /„ r .a, 

%\U ae A £ifl B(a,»?) 

In view of Lemma 17.51 letting rj — > 0, we obtain 



p£ [u e ] < ^ (2 / iV^l^x' + vr^lA^llnrz + ^IA^I ) + o £ (l) + o v (l) 



| W(^) dPfc) + O e (l) + Ofl(l)) + ^1^1 + O, 



F e [u e ] < f 4 / W(s>) dP(<p) + o e (l) + 0^(1) + ^^|A £ ,^| + o e (l). 



Letting e — > 0, using that r £ — )• 3 1 / 3 and the fact that |A £) h| = ^m Ej R\£ £ \ 2 with m et R 
and then finally letting R — )■ oo, we conclude that 



m, 



limsu P F £ [n e ] < 3 4/3 / W^dP^) + 

J 



3 4/3 



J?? 



Since |3 2 / 3 m = |(<5 — 8 C ), this completes the proof of part ii) of Theorem [TJ □ 
7.4 Proof of Theorem [2] 

In order to prove Theorem [2J it suffices to show that 

min F°[P] = 3 4 / 3 min W + - ^ . (7.23) 

For the proof, we use the following result, adapted from Corollary 4.4 in [35] . 

Proposition 7.7 (Corollary 4.4 in |35j). Let <p € „4,i 6e given, such that W((p) < oo. For 
any R such that R 2 £ 27rN, there exists a 2R-periodic ipR such that 

-Aipn = 2vr EaGA fl S a ~ 1 in 

— — = on dK R , 

ou 

where Ar is a finite subset of the interior of Kr, and such that 

W(<pr,1k r ) . „ rl v 
hmsup —— < W(<p). 

Let us take ip to be a minimizer of W over A m (which exists from [35]). We may 
rescale it to be an element of Ai- Then Proposition 17. 71 yields a (fR, which can be extended 
periodically. We can then repeat the same construction as in the beginning of this section, 
starting from VipR instead of bR, and in the end it yields a u £ with 

lim sup F £ [u e ] < 3 4/3 min W + t!l^_lM. . 

e _>0 ■Am 8 
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It follows that 

lim sup min F £ < 3 4 / 3 min W H — . 

E— >0 A m 8 

But by part i) of Theorem [1] applied to a sequence of minimizers of F £ , we also have 

lim inf min F £ > inf F° > 3 4/3 min W + 3 ' ^ - ^ 
e^O A V A m 8 

where the last inequality is an immediate consequence of the definition of F°. Comparing 
the inequalities yields that there must be equality and (17.23|) is proved, which completes 
the proof of Theorem [2J □ 
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